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Abstract. An overview of the author's results is given. Property C stands 
for completeness of the set of products of solutions to homogeneous linear 
Sturm-Liouville equations. The inverse problems discussed include the clas- 
sical ones (inverse scattering on a half-line, on the full line, inverse spectral 
problem), inverse scattering problems with incomplete data, for example, in- 
verse scattering on the full line when the reflection coefficient is known but no 
information about bound states and norming constants is available, but it is a 
priori known that the potential vanishes for x < 0, or inverse scattering on a 
half-line when the phase shift of the s-wave is known for all energies, no bound 
states and norming constants are known, but the potential is a priori known 
to be compactly supported. If the potential is compactly supported, then it 
can be uniquely recovered from the knowledge of the Jost function f{k) only, 
or from /'(O, fc), for all fc £ A, where A is an arbitrary subset of (0,oo) of 
positive Lebesguc measure. 

Inverse scattering problem for an inhomogeneous Schrodinger equation is 
studied. 

Inverse scattering problem with fixed-energy phase shifts as the data is 
studied. 

Some inverse problems for parabolic and hyperbolic equations are investi- 
gated. 

A detailed analysis of the invertibility of all the steps in the inversion pro- 
cedures for solving the inverse scattering and spectral problems is presented. 

An analysis of the Newton-Sabatier procedure for inversion of fixed-energy 
phase shifts is given. 

Inverse problems with mixed data are investigated. 

Representation formula for the /-function is given and properties of this 
function are studied. 

Algorithms for finding the scattering data from the /-function, the /- 
function from the scattering data and the potential from the /-function are 
given. 

A characterization of the Weyl solution and a formula for this solution in 
terms of Green's function are obtained. 
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1. Property C for ODE 

In this paper a review of the author's results is given and some new results are 
included. The bibliography is not complete. Only the papers and books used in 
the work on this paper are mentioned. The contents of this paper are clear from 
the table of contents. 

The results presented in this paper include: 

1. Property C for ordinary differential equations (ODE), that is, theorems about 
completeness of the sets of products of solutions to homogeneous ODE. 

2. Uniqueness theorems for finding the potential 

a) from the /-function (which equals the Weyl function), 

b) from the classical scattering data for the half-axis problem (a new very 
short proof which does not use the Marchenko method) , 

c) from the phase shift of s-wave in the case when the potential q is com- 
pactly supported and no bound states or norming constants are known, 

d) from the reflection coefficient only (when q = for x < xq), 

e) from mixed data: part of the set of eigenvalues and the knowledge of 
q{x) on a part of a finite interval, 

f) from overdetermined Cauchy data, 

g) from part of the fixed-energy phase shifts, 

h) from various type of data which are typical in PDE problems, 

i) from f(k) or /'(O, k) when q is compactly supported, 

j) from the scattering data for a solution to an inhomogeneous Schrodinger 
equation. 

3. Reconstruction algorithms for finding the potential from overdetermined Cauchy 
data, for finding f{k) and f'{0,k) from the scattering data, for finding the 
scattering data from the /-function and the /-function from the scattering 
data. 

4. Properties of the /-function and a representation formula for it. 

5. Stability estimate for the solution of the inverse scattering problem with fixed- 
energy data. Example of two compactly supported real-valued piecewise- 
constant potentials which produce practically the same phase shifts for all 
values of £. 

6. Discussion of the Newton-Sabatier procedure for inversion of the fixed-energy 
phase shifts. Proof of the fact that this procedure cannot recover generic 
potentials, for example, compactly supported potentials. 

7. Detailed analysis and proof of the invertibility of each of the steps in the 
inversion schemes of Marchenko and Gelfand-Levitan. 

8. Representation of the Weyl solution via the Green function and a characteri- 
zation of this solution by its behavior for large complex values of the spectral 
parameter and x running through a compact set. 

Completeness of the set of products of solutions to ODE has been used for inverse 
problems on a finite interval in the works of Borg j|] and Levitan , . 

Completeness of the set of products of solutions to homogeneous partial dif- 
ferential equations (PDE) was introduced and used extensively under the name 
property C in |^-[^, and Property C in these works differs essentially from 
the property C defined and used in this paper: while in [p6[-[^, and JTzt property 
C means completeness of the set of products of solutions to homogeneous PDE with 
fixed value of the spectral parameter, in this paper we prove and use completeness 
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of the sets of products of solutions to homogeneous ordinary differential equations 
(ODE) with variable values of the spectral parameter. Note that the dimension of 
the null-space of a homogeneous PDE (without boundary conditions) with a fixed 
value of the spectral parameter is infinite, while the dimension of the null-space 
of a homogeneous ODE (without boundary conditions) with a fixed value of the 
spectral parameter is finite. Therefore one cannot have property C for ODE in the 



sense of [^-|3l|, and |g/7|, because the set of products of solutions to homogeneous 
ODE with fixed value of the spectral parameter is finite-dimensional. 

In this paper property C for ordinary differential equations is defined, proved 
and used extensively. Earlier papers are Q and Q . 

Let 

(I.l) (u := u" + k'^u - q{x)u = 0, a:: G M = (-cx), oo). 

Assume 



(1.2) qeLis, Li, 



/OO 
il + \x\r\q{x)\dx<^}. 
-OO 



It is known Q, ||T^ that there is a unique solution (the Jost solution) to (I.l) 
with the asymptotics 

(1.3) /(a;,fc) = e'*=^ + o(I), x +oo. 

We denote f+{x,k) := f{x,k), f-{x,k) f{x,-k), fc e M. The function /(O, fc) = 
f{k) is called the Jost function. The function f{k) is analytic in C+ := {fc : Imk > 
0} and has at most finitely many zeros in C+ which are located at the points 
ikj,kj > 0, 1 < j < J. The numbers — fc| are the negative eigenvalues of the 

selfadjoint operator defined by the differential expression Lg := —-^ -\- q{x) and 
the boundary condition u{0) = in L^(K+), R+ = (0,00). The function /(fc) may 
have zero at fc = 0. This zero is simple: if /(O) = then /(O) 7^ 0, / := 
Let ip and ip be the solutions to (I.l) defined by the conditions 

(1.4) ^(0,fc) = 0, <y9'(0,fc) = l; V'(0,fc) = l, ?/''(0,fc) = 0, 

where ip' :— It is known |Q, |jl^, that ip{x,k) and ip{x,k) are even entire 
functions of fc of exponential type < 

Let g±{x,k) be the unique solution to (1.1) with the asymptotics 

(1.5) 5±(a;,fc) = e*^'^^ + 0(1), a; ^ -00 
Definition 1.1. Let p{x) E Li,i(]R+) and assume 

(1.6) / p{x)h{x,k)j2{x,k)dx = Q, Vfc>0, 



Jo 

where fj{x,k) is the Jost solution to (1.1) with q{x) ~ qj{x), j ~ 1,2. // (1.6) 

implies p{x) — 0, then we say that the pair {Li, Lj := Lq^ := '^Qji^) 
property C+. 

Ifp € ii4(R-) and 

(1.7) / p{x)gi{x,k)g2{x,k)dx = Q Vfc > 0, 

implies p{x) — 0, then we say that the pair {Li, L2} has property C-. 
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In (1.7) gj gj+{x,k). 

Fix an arbitrary b > 0. Assume that 

(1.8) / p{x)ipi{x,k)(p2ix,k)dx ^0 Vfc>0 
Jo 

implies p{x) = 0. Then we say that the pair L2} has property and similarly 
is defined, ipj replacing ipj in (1.8). 

Theorem 1.1. The pair {Li,L2\, where Lj := —-^ +qj{x), Qj £ Li i(M+),j = 
1,2, has properties C+,C^p and C^. If qj G Li_i(R_), then {ii,L2} has property 

Proof. Proof can be found in [|l8j . We sketch only the idea of the proof of property 
Using the known formula 

(1.9) f,{x,k) = e'''- + Aj{x,y)e'''ydy, j = l,2, 



where Aj{x,y) is the transformation kernel corresponding to the potential qj{x), 
j = 1,2, see also formula (2.17) below, and substituting (1.9) into (1.6), one gets 
after a change of order of integration a homogeneous Volterra integral equation for 
p{x). Thus p{x) =0. □ 

The reason for taking & < 00 in (1.8) is: when one uses the formula 

/ ,x s'm(kx) , ,sin{ky) , ^ „ 

(1.10) ^^.(a:,/c) = — ^+y^ if,(x,y)— ^dy, j = l,2, 

for the solution Lpj to (1.1) (with q — qj) satisfying first two conditions (1.4), 
then the Volterra-type integral equation for p{x) contains integrals over the infinite 
interval {x,oo). In this case the conclusion p{x) = does not hold, in general. 
If, however, the integrals are over a finite interval {x,b), then one can conclude 
p{x) = 0. 

The same argument holds when one proves property C^, but formula (1.10) is 
replaced by 



(1.11) ^j{x,k) = cos{kx) + / Kj{x,y) cos{ky)dy, j = 1,2, 

Jo 

with a different kernel Kj(x,y). 

2. Applications of property C 

2.1. Uniqueness of the solution inverse scattering problem with the data 
I-function. The I-function I{k) is defined by the formula 

This function is equal to the Weyl function m{k) which is defined as the function 
for which 

(2.2) W{x,k) ■.= ^p{x,k)+m{k)Lp{x,k) e L^{M.+ ), Imk > 0, 
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where W{x,k) is the Weyl solution, W{0,k) = 1, W'{0,k) = m{k). Note that 
W{x,k) — as follows from formulas (1.3), (2.1) and from formula (2.3) 

which says I(k) = m(k). 
To prove that 

(2.3) I{k)=m{k), 

one argues as follows. If g G Li_i(R+), then there is exactly one, up to a constant 
factor solution to (1.1) belonging to L^(R+) when Imk > 0. 
Since f{x, k) is such a solution, one concludes that 

(2.4) f{x,k) = c{k)[^P{x,k) +m{k)ip{x,k)], c{k) ^ 0. 
Therefore, 

_ f{0,k) + mik)^'{0,k) 
^ '^ mk)+m{k)i;{0,k) ^ ^ 

as claimed. 

In sections 11 and 12 the /-function is studied in more detail. 
The inverse problem (IPl) is: 
Given I{k) for all fc > 0, find q{x). 

Theorem 2.1. The IPl has at most one solution. 



Proof. Theorem 2.1 can be proved in several ways. One way [[l9[ is to recover the 
spectral function p(A) from /(fc), k — . This is possible since Iml{k) — ^ 
k > 0, and 



(2-5) dp{X} = {:}fj^^)\^' 



A > 0, 



Ej=i Cj'5(A + fc2) dA, A < 0, kj > 0, 

where —k^ are the bound states of the Dirichlet operator Lq = + q{x) in 

L'^(R+), f{ikj) — 0, (5(A) is the delta-function, and 

Note that ikj and the number J in (2.5) can be found as the simple poles of I{k) 
in C+ and the number of these poles, and 

(2.7) Cj = -2ikj Re^s I{k) = 2kjrj, 

where ir^ := ReSk=kj I{k), so = 

It is well known that dp{X) determines q{x) uniquely ll4[], [^ . An algorithm for 

recovery of q{x) from dp is known (Gelfand-Levitan) . In [|l9| a characterization of 

the class of /-functions corresponding to potentials in C™^ (R-i- ) , m > is given. 
Here we give a very simple new proof of Theorem 2.1 (cf flql ): 
Assume that qi and q^ generate the same /(fc), that is, Ii{k) = hik) := I{k). 

Subtract from equation (1.1) for fi(x,k) this equation for f2{x,k) and get: 

(2.8) Liw=pf2, p(x) := qi(x) - q2(x), w -.^ fi{x, k) - f2{x, k). 
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Multiply (2.8) by /i and integrate by parts: 

p{x)f2{x,k)h{x,k)dx^{w'h~wf[)\^ = (/l/2-/{/2)L=o 

(2.9) = hUhik) - hm = VA; > 0, 

where we have used (1.3) to conclude that at infinity the boundary term vanishes. 
From (2.9) and property C+ (Theorem 1.1) it follows that p{x) — 0. Theorem 2.1 
is proved. 

□ 

2.2. Uniqueness of the solution to inverse scattering problem on the half 
axis. This is a classical problem jljl, [0. The scattering data are 

(2.10) S = {S[k), kj, s„ l<j< J}. 
Here 

(2.11) Sik) 

is the S-matrix, kj > are the same as in section 2.1, and the norming constants 
Sj are the numbers 

(2.12) s, := — : > 0. 

/(*fc,)/'(0,z%) 

Note that (2.6) implies 

(2.13) cj=s,[f{0,tkj)]\ 

Theorem 2.2. Data (2.10) determine q{x) e ii_i(M+) uniquely. 

Proof. This result is due to Marchenko |l4| . We give a new short proof based on 
property C (fl^j). We prove that data (2.10) determine I{k) uniquely, and then 
Theorem 2.2 follows from Theorem 2.1. To determine I{k) we determine f{k) and 
/'(0,fc) from data (2.10). 

First, let us prove that data (2.10) determine uniquely f{k). Suppose there are 
two different functions f{k) and h{k) with the same data (2.10). Then 

("^' m-wwy 

The left-hand side in (2.14) is analytic in since f{k) and h{k) are, and the 
zeros of h{k) in C+ are the same as these of /(fc), namely «fcj, and they are simple. 
The right-hand side of (2.14) has similar properties in C_. Thus is an entire 
function which tends to 1 as |fc| oo, so, — 1 and /(fc) = h{k). The relation 

(2.15) lim f{k) = 1 

foUows from the representation 

(2.16) /(fc) = 1 + / A(0, y)e'''ydy, A(0, y) e Li(R+). 
Various estimates for the kernel A{x, y) in the formula 

(2.17) f{x,k) = e^'^- + A{x,y)e''^ydy 
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are given m 
(2.18) 

(2.19) 

(2.20) 



We mention the following: 
\A{x,y)\<ca > ^ ^ 



dA{x,y) 1 / x + y 
dx 4'^^^ 2 

dA{x,y) 1 f x + y 



a{x) := 

< ca{x)a 

< ca{x)a 



\qit)\dt, 

x + y '' 
2 

x + y 



where c > here and below stands for various estimation constants. 
From (2.17) and (2.18) formula (2.16) follows. 
Thus, we have proved 

(2.21) S^fik). 
Let us prove 

(2.22) S^f'(Q,k). 
We use the Wronskian: 

(2.23) no, k)f{-k) - /'(O, ~k)f{k) = 2ik, fc e M. 



The function /(fc) and therefore f{—k) = f{k), where the overbar stands for 
complex conjugate, we have already uniquely determined from data (2.10). Assume 
there are two functions /'(O, k) and h'{0, k) corresponding to the same data (2.10). 
Let 

(2.24) w{k) -.^ f{Q,k)- h'{0,k). 

Subtract (2.23) with h'{0, ±fc) in place of /'(O, ±fc) from equation (2.23) and get 

w{k)f{-k)~w{~k)f{k) = 0, 

or 

w{k) w{—k) 



(2.25) 



f{k) f{-k) 



Vfc e 



Claim, jj^ is analytic in C+ and vanishes at infinity and is analytic in 

C- and vanishes at infinity. 



"(-fc) 



If this claim holds, then 



0, fc € C, and therefore w{k) = 0, so /'(O, k) = 

h'{0,k). 

To complete the proof, let us prove the claim. 
From (2.17) one gets: 



(2.26) 

Taking k 

(2.27) 

Thus 
(2.28) 



/'(O, k)^ik- A{Q, 0) + / A,(0, y)e''^ydy 



-oo in (2.16), integrating by parts and using (2.20), one gets: 



ik 



1 f 

-- Ay{Q,k)e'''y dy. 
ik Jo 



A{Q,0) = - lim ik[f{k) - 1]. 

k — >oo 
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Since f{k) is uniquely determined by data (2.10), so is the constant A(0,0) (by 
formula (2.28)). 

Therefore (2.24) and (2.28) imply: 

(2.29) lim w(k) = 0. 

\k\^oo, keC+ 

It remains to be checked that (2.10) implies 

(2.30) w{ikj) = 0. 

This follows from formula (2.12): if f{k) and sj are the same, so are f'{0,ikj), and 
w{ikj) = as the difference of equal numbers 

l-ik ■ 

h'{0,ikj) = f'{0,ikj) = — ' 



f{ikj)sj 

Theorem 2.2 is proved. □ 

In this section we have proved that the scattering data (2.10) determines the I- 
function (2.1) uniquely. The converse is also true: implicitly it follows from the 
fact that both sets of data (2.1) and (2.10) determine uniquely the potential and 
are determined by the potential uniquely. A direct proof is given in section 12 
below. 

2.3. Compactly supported potential is uniquely determined by the phase 
shift of s-wave. Consider the inverse scattering on half-line and assume q{x) = 
for a; > a > 0, where a > is an arbitrary fixed number. 

The phase shift of s-wave is denoted by 6{k) and is defined by the formula 

(2.31) f{k) = |/(fc)|e-^^W, 
so the S-matrix can be written as 

(2.32) 5(fc) = 4^=e2'*W. 

/(«) 

If q{x) is real-valued, then 

(2.33) 5{-k) = -S{k), keR, 
and if g e Li,i(]R+), then 

(2.34) 5{oo) = 0. 
Note that S-matrix is unitary: 



(2.35) S{-k) = S{k), \S{k)\ = 1 if A; e M. 
Define index of S{k): 

(2.36) v — indSik)-- — I d\nS{k) = —A^arg S{k). 
From (2.32), (2.33) and (2.34) one derives a formula for the index: 

y =-An5{k) = -[(5(-0) - ^(-oo) + 5{+oo) - 5{+Q)] 

TT TT 

(2-37) _2 r-2J if/(0)^0, 

TT ^ ' [-2J-1 if/(0)=0. 
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Here we have used the formula: 

(2.38) -S{+Q) = #{zeros of ,f{k) in C+} + i<5o, 

TT 2 

which is the argument principle. Here So :— if /(O) ^ and Sq 1 if /(O) = 0. 
The zero of /(fc) at fc = is called a resonance at zero energy. 
Let us prove the following result [ p^ : 

Theorem 2.3. Ifq £ Li i{R^) decays faster than any exponential: \q{x)\ < ce~'^'^'^, 
7 > 1, then the data {S{k) Vfc > 0} determines q{x) uniquely. 

Proof. Our proof is new and short. We prove that, if q is compactly supported 
or decays faster than any exponential, e.g. \q{x)\ < ce"^'^'^, 7 > 1, then S{k) 
determines uniquely kj and Sj, and, by Theorem 2.2, q{x) is uniquely determined. 

We give the proof for compactly supported potentials. The proof for the poten- 
tials decaying faster than any exponentials is exactly the same. The crucial point 
is; under both assumptions the Jost function is an entire function of k. 

li q{x) is compactly supported, q{x) — for x > a, then /(fc) is an entire function 
of exponential type < 2a, that is |/(fc)| < ce^'^l'^'l p. 278]). Therefore 5(fc) 

is meromorphic in C+ (see (2.32)). Therefore the numbers kj, 1 < j < J, can be 
uniquely determined as the only poles of S{k) in C+. One should check that 

(2.39) f{-tkj)^0 if /(*fc,) =0. 

This follows from (2.23): if one takes k = ikj and uses f{ikj) ~ 0, then (2.23) 
yields 

(2.40) /'(O, ikj)J{~ikj) = -2kj < 0. 

Thus f{—ikj) 0. Therefore d{k) determines uniquely the numbers kj and J. 
To determine Sj, note that 

(2.41) Res 5(fc) = /HM ^ i,^, 

as follows from (2.12) and (2.40). Thus the data (2.10) are uniquely determined 
from S{k) if q is compactly supported, and Theorem 2.2 implies Theorem 2.3. □ 

Corollary 2.1. If q{r) G iyi^i(M+) is compactly supported, then the knowledge of 
f{k) on an arbitrary small open subset 0/ M+ (or even on an infinite sequence 
kn > 0, fc„ 7^ km if m n, kn k as n 00, k > 0) determines q{r) uniquely. 
In section 4 we prove a similar result with the data f'{0,k) in place of f{k). 

2.4. Recovery of q E Li i(R) from the reflection coeflicient alone. Consider 
the scattering problem on the full line: u solves (1.1) and 

(2.42) u ^ t{k)e''''= , x^+00, 

(2.43) u-e*'="+r(/fc)e-''=", x ^ -00. 

The coefficients t{k) and r{k) are called the transmission and reflection coeffi- 
cients (see |l^] and [0). In general r{k) alone cannot determine q{x) uniquely. 
We assume 

(2.44) q{x) = for a; < 0, 

and give a short proof, based on property C, of the following: 
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Theorem 2.4. If q £ Li,i(M) and (2.44) holds, then r{k), Vfc > 0, determines 
q{x) uniquely. 

Proof. We claim that r{k) determines uniquely I{k) if (2.44) holds. Thus, Theorem 
2.4 follows from Theorem 2.1. To check the claim, note that u{x,k) = t{k)f{x,k), 
so 

and use (2.44), (2.43) to get u = e**"^ + r{k)e^^^'' for x < 0, so 

u'{Q,k) ^ ik{l-r{k)) 
^ ■ ' u{0,k) l + r{k) ■ 

Prom (2.45) and (2.46) the claim follows. Theorem 2.4 is proved. □ 

2.5. Inverse scattering with various data. Consider scattering on the full line 
(1.1), (2.42)-(2.43), assume q{x) = if a; [0, 1], and take as the scattering data 
the function 

(2.47) u{0,k) :=uo{k), Vfc > 0. 

Theorem 2.5. Data (2.47) determine q{x) uniquely. 

Proof. If (7 = for a; < 0, then u{x, k) = e"^"^ + r(fc)e~'''^ for x < 0, m(0, k) = 
1 + r(k), so data (2.47) determine r(k) and, by Theorem 2.4, q{x) is uniquely 
determined. Theorem 2.5 is proved. Of course, this theorem is a particular case of 
Theorem 2.4. □ 



Remark 2.1. Other data can be considered, for example, u'{0,k) := v{k). Then 
u'{0,k) = ik[l — r{k)], and again v{k) determines r{k) and, by Theorem 2.4, q{x) 
is uniquely determined. 

However, if the data arc given at the right end of the support of the potential, 
the inverse problem is more difficult. For example, if u(l,k) := ui{k) is given for 
all A; > 0, then ui{k) = t{k)e'^'', so t{k) is determined by the data uniquely. 

The problem of determining q{x) from t{k) does not seem to have been studied. 
If q{x) > 0, then the selfadjoint operator Lg = — ^ +q{x) in L'^{R) does not have 
bound states (negative eigenvalues). 

In this case the relation |r^(A;)| + |t^(fc)| = 1 allows one to find 

\r{k)\ = Vl-|t(fc)P, kGR. 

Define 

(2.48) aik):=expl-^ r'^ds}. 

I T^lJ-oo S-k } 

The function (2.48) has no zeros in C+ if Lg has no bound states. If q is compactly 
supported then r{k) and t{k) are meromorphic in C. 
Let us note that 

(2.49) f{x,k) := f+{x,k) =b{k)g-{x,k) + a{k)g+{x,k) 



(2.50) 



g-{x,k) = c{k)f+{x,k) + d{k)f_{x,k). 
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It is known 0, that 

(2.51) a(-fc)=^I(fc), b{~k)=b(k), keR, 

(2.52) c{k) = -b{-k), d{k) = aik), 

(2.53) a{k)^-^[f{x,k),g^(x,k% b{k) = ^[f+{x,k), g+{x,k)], 

where [f,g] '■— fg' — f'g is the Wronskian, 

(2.54) \a{k)\^ = l + \b{k)\\ 

The function a{k) is analytic in C+. One can prove |l4| p. 288] 

r°° q(x)dx /l\ 

(2.56) «(fc) = 1 - 
and 

(2.57) 5(fc) = o(^^^, /j^c». 

The function r(fc) does not allow, in general, an analytic continuation from the 
real axis into the complex plane. However, if q{x) = for a; < 0, then b{k) admits 
an analytic continuation from the real axis into C+ and r{k) is meromorhic in C+. 

If q{x) is compactly supported the functions f±{x,k) and g±{x,k) are entire 
functions of k of exponential type, so that r(fc) and tik) are meromorphic in C. 
From (2.54) one can find \b{k)\ since a{k) is found from formula (2.48) (assuming 
no bound states). 

The conclusion is: recovery of a compactly supported potential from the trans- 
mission coefficient is an open problem. 

3. Inverse problems on a finite interval with mixed data 

Consider equation (1.1) on the interval [0,1]. Take some selfadjoint boundary 
conditions, for example: 

(3.1) ~u" + q{x)u-\u^Q, 0<a;<l; u(0) = m(1) = 0, A = fc^ 

Assume q £ L^[0, l],g = q. Fix b € (0,1) arbitrary. Suppose q{x) is known on 
the interval [6, 1] and the subset {Xm(n)} of the eigenvalues of the problem (3.1) is 
known, n — 1,2,... where m{n) is a sequence with the property 

(3.2) !!!M = i(i+e„), e„^0, a > 0. 

n a 

Theorem 3.1. The data {Am(n),'T' = I7 2, (7(x), 6 < a; < 1} determine uniquely 
9(2^), < X < b, if a > 2b. If a — 2b and J2'^=i kl ^ then the above data 
determine q{x), < x < b, uniquely. 
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Proof. First, assume a > 2b. If there are qi and 52 which produce the same data, 
then as above, one gets 

(3.3) 



G{X) := g{k) := / p{x)ipi{x, k)(p2{x, k) dx = {ipiw' — ip'iw) = {tfiw' — tp'^w) 
Jo 

where w :— ipi ~ ip2, p := qi — q2, k = VA- Thus 

(3.4) g{k) = at fc = ±^ X,n{n) ±fcn. 

The function G(A) is an entire function of A of order i (see formula (1.10) with 
k — \/A), and is an entire even function of k of exponential type < 2b. One has 

2h\Imk\ 

(3.5) I^WI^^TW 
The indicator of g is defined by the formula 

(3.6) h{0) := hJe) TSr^^^^^^^, 

where k = re*^. Since |/m/c| = r\ sin0|, one gets from (3.5) and (3.6) the following 
estimate 

(3.7) /i(6l) < 26|sin6l|. 

It is known jlO| formula (4.16)] that for any entire function g{k) ^ of expo- 
nential type one has: 

n(r) 1 f^" 

(3.8) liniAi<— / hg{9)d9, 

where n(r) is the number of zeros of g(k) in the disk \k\ < r. From (3.7) one gets 

1 f"^^ 2b f^^ 46 

(3.9) — / hg{e)d9<— \sm9\de = - 



From (3.2) and the known asymptotics of the Dirichlet eigenvalues: 

(3.10) A.,i = (7m)^ + c + 0(1), n — > 00, c — const, 
one gets for the number of zeros the estimate 

(3.11) n(r)>2 ^ 1 = 2 — [1 + o(l)], r 00. 

From (3.8), (3.9) and (3.11) it follows that 

(3.12) a < 2b. 

Therefore, if a > 26, then g{k) = 0. If g{k) = then, by property dp (Theorem 
1.1), p{x) — 0. Theorem 3.1 is proved in the case a > 2b. 
Assume now that a = 2b and 

00 

(3.13) ^|£„|<c3o. 

We claim that if an entire function G(A) in (3.3) of order ^ vanishes at the points 

2 ? 

(3.14) il + Sn), 
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and (3.13) holds, then G(A) = 0. If this is proved, then Theorem 3.1 is proved as 
above. 

Let us prove the claim. Define 



(3.15) 

and recall that 
(3.16) $o(A) 



$(A) 
sin((T\/A) 



A 



Since G(A„) = 0, the function 



n—l ^ 

oc ^ 

n 1 



A 



(3.17) 



u;(A) 



GjX) 



is entire, of order < i. Let us use a Phragmen-Lindelof lemma. 

Lemma 3.1. pO| , Theorem 1.22] // an entire function w{X) of order < 1 has the 
property sup^oo<y<oo\w{iy)\ < c, then w{X) = c. If, in addition w(iy) — s- as 
y — > +00, then w{X) = 0. 

We use this lemma to prove that w{X) = 0. If this is proved then G(A) = and 
Theorem 3.1 proved. 

The function w(A) is entire of order ^ < 1. 
Let us check that 



(3.18) 



sup \w[iy)\ < oo, 

-oo<y<oo 



and that 

(3.19) \w{iy)\ — > as y ^ +oo. 

One has, using (3.5), (3.15), (3.16) and taking into account that a = 2b: 
(3.20) 



G{iy) $o(«y) 



$(iy) $o(iy) 



< 



^2b\Im-/vy\ I ^a\Imy/%y\ 



(l + |y|) 



1 



\y\ 



n , , 

h=i \i 



<- 



n 



Xn 



< 



1 



\y\ 



n (1 



kJ) < 



Cl 



Here we have used elementary inequalities: 
(3.21) 



1 + a a „ , 1 + a 

< - if a > d > 0; : < 1 if < a < d, 

1 + d ~ d ~ 1 + d ~ 



with a := d := -p-, and the assumption (3.13). 



We also used the relation: 

sm{a\/Jy) 



as y 



'OO. 



Estimate (3.20) implies (3.18) and (3.19). An estimate similar to (3.20) has been 
used in the literature (see e.g.[Q). 

Theorem 3.1 is proved. □ 
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Remark 3.1. Theorem 3.1 yields several results obtained in [Q, and an earlier 
result of Hochstadt-Lieberman which says that the knowledge of all the Dirichlet 
eigenvalues and the knowledge of q{x) on [i, l] determine uniquely q{x) on [O, i] . 

In this case b = i, a = 1. 

One can also obtain the classical result of Borg Q and its generalization due to 
Marchenko 

Two spectra (with the same boundary conditions on one of the ends of the 
interval and two different boundary conditions on the other end) determine q{x) 
and the boundary conditions uniquely. 

4. Property C and inverse problems for some PDE 
4.1. Consider the problem 

(4.1) ut = - q{x)u, < a; < 1, <> 0, 

(4.2) M(a;,0) = 0, 

(4.3) u(0,t) = 0, u{l,t)=a{t). 

Assume the a{t) ^ is compactly supported, a{t) S L^(0,oo), q{x) £ i^[0, 1], 
problem (4.1) - (4.3) is solvable, and one can measure the data 

(4.4) u'(l, t) u^{l, t) := b{t) Mt > 0. 

The inverse problem (IPl) is: 
Given {a(i), b{t), 'it > 0} find q{x). 

Theorem 4.1. IPl has at most one solution. 

Proof. Laplace-transform (4.1) - (4.3) to get 

j>OC 

(4.5) v" - \v - q{x)v ^ Q < a; < 1, w := / e^^*u{x,t) dt, 

Jo 

(4.6) v{0,X)=0, v{l,X) ^ A{X) := e-^^a{t)dt. 



(4.7) v'{l,\) = B{\):= b{t)e-^'dt. 

Jo 

Assume that there are qi{x) and q2{x) which generate the same data {A{X), B{X), 
VA > 0}. hct p{x) := qi{x) — q2{x),w := VI—V2. Subtract from equation (4.5) with 
w = Vi, q = qi, similar equations with v = V2, q ~ q2, and get 

(4.8) iiw :— w" — Xw — qiw — pv2, 

(4.9) w(0,A)=0, u;(l, A) = w'(l, A) 0. 

Multiply (4.8) by ipi{x, A), where iiipi = 0, i^i(0. A) = 0, A) = 1, integrate 

over [0, 1] and then by parts on the left-hand side, using (4.9). The result is: 

(4.10) / p{x)v2{x,X)Lpi{x,X)dx^Q VA > 0. 



Note that '^i{x, A) is an entire function of A. 
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Since a{t) ^ and is compactly supported, the function A{X) is an entire function 
of A, so it has a discrete set of zeros. Therefore V2(x,X) = c{X)lp2{x; X) where 
c(A) ^ for almost all A e M+. 

Property (Theorem 1.1) and (4.10) imply p{x) = 0. Theorem 4.1 is proved. 

□ 

Remark 4.1. One can consider different self adjoint homogeneous boundary condi- 
tions at X = 0, for example, m'(0, t) ~ or u'(0, t) — /iou(0, t) = 0, ho — const > 0. 

A different method of proof of a result similar to Theorem 4.1 can be found in 
in and in g. In |] some extra assumptions are imposed on q{x) and a(t). 



4.2. Consider the problem: 

(4.11) utt ^ Uxx - q{x)u, a; > 0, t > 0, 

(4.12) M = ut = at i = 0, 

(4.13) u(0,t) = (5(t), 
where 5{t) is the delta-function. Assume that 

(4.14) q{x) = Q if x>l, q^q, q(iL^[0,l]. 
Suppose the data 

(4.15) u{l,t):=a{t) 

are given for alH > 0. 

The inverse problem (IP2) is: 
Given a{t) Vt > 0, find q(x). 

Theorem 4.2. The IP2 has at most one solution. 

Proof. Fourier-transform (4.11)-(4.13), (4.15) to get 

(4.16) v" + k'^v - q{x)v x>Q, 

(4.17) «(0,fc) = l, 

(4.18) v(l,k)^A{k)~ a{t)e'^'dt, 

Jo 

where 

(4.19) v{x,k) = / u{x,t)e''''dt. 

Jq 

It follows from (4.19) and (4.16) that 

(4.20) v{x,k)^cik)fix,k), 

where f{x,k) is the Jost solution to (4.16). From (4.20) and (4.17) one gets 

(4.21) "(-^-fr- 

where f{k) = /(O, fc). From (4.21) and (4.18) one obtains 
(4.22, fik) ^ 41^. 
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Prom (4.14) and (4.16) one concludes 

(4.23) /(a;,A;)=e''=^ fora;> 1. 
Prom (4.23) and (4.22) one gets 

(4.24) m = 

Thus f{k) is known for all k > 0. 

Since q{x) is compactly supported, the data {f{k), \/k = 0} determine q{x) 
uniquely by Theorem 2.3. Theorem 4.2 is proved. 

□ 



RemEirk 4.2. One can consider other data at x = 1, for example, the data Ux{l, t) 
or Ua;(l,t) + hu{l,t). The argument remains essentially the same. 

However, the argument needs a modification if (4.13) is replaced by another 
condition, for example, Ux{0^t) = 5{t). In this case v{x,k) = jr^^, and in place 
of f{k) one obtains /'(0,fc) Vfc > from the data (4.18). 

The problem of finding a compactly supported q{x) from the data {/'(O, k) VA; > 
0} was not studied, to our knowledge. We state the following: 

Claim. The data /'(O, fc) known on an arbitrary small open subset of (0, oo) or 
even on an infinite sequence of distinct positive numbers A;„ which has a limit point 
k > 0, determines a compactly supported q{r) G L^(M+) uniquely. 

Our approach to this problem is based on formula (2.23). If f'{0,k) is known 
for all fc > 0, then /'(0,-fc) = /'(0,fc) is known for all fc > 0, and (2.23) can be 
considered as the Riemann problem for finding /(fc) and /(— fc) from (2.23) with the 
coefficients /'(O, fc) and /'(O, — fc) known. If ^(a;) G Li(M+) is compactly supported 
then /'(O, fc) is an entire function of fc. Thus the data determine /'(O, fc), for all 
fc > 0. 

We want to prove that (2.23) defines /(fc) uniquely if /'(O, fc) is known for all 
fc > 0. Assume the contrary. Let /(fc) and h{k) be two solutions to (2.23), and 
w := f — h, w{k) — » as |fc| ^ oo, fc e C+. Then (2.23) implies 

The function f'{0,k) has at most finitely many zeros in C+. All these zeros 
are at the points iKj, 1 < j < Ji, where — k| arc the negative eigenvalues of the 

Neumann operator Lg := --£^ + q{x) in L2(R+), u'(0) = 0. Also /'(0, 0) may 
vanish. 

From (2.23) one concludes that w{iKj) = if f'(0,iKj) = 0. Indeed, one has 
w{iKj)f' {0, —iKj) = iv(—iKj)f'{0,iKj). If f'{Q,iKj) = then /'(O, — iKj) 7^ as 
follows from (2.23). Therefore w{iKj) = as claimed, and the function ^Ti^qI-j is 
analytic in C-|- and vanishes at infinity in C+. Similary, the right-hand side of (4.25) 
is analytic in C_ and vanishes at infinity in C_. Thus, by analytic continuation, 

u(k) 



is an entire function which vanishes at infinity and therefore vanishes iden- 
tically. Therefore w{k) = and /(fc) = h{k). Thus, the data {/'(0,fc), Vfc > 0} 
determines uniquely {/(fc), Vfc > 0}. 
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Since q{x) is compactly supported, Theorem 2.3 implies that q{x) is uniquely 
determined by the above data. The claim is proved. □ 

5. Invertibility of the steps in the inversion provedures in the 
inverse scattering and spectral problems 

5.1. Inverse spectral problem. Consider a sclfadjoint operator Lg in _L^(R+) 
generated by the differential expression Lq = + q(x), q{x) E Lj^^(R.+ ), q{x) = 
q{x)^ and a selfadjoint boundary condition at x = 0, for example, m(0) = 0. Other 
selfadjoint conditions can be assumed. For instance: u'(0) = /low(O), ho = const > 
0. 

We assume that q{x) is such that the equation (1.1) with Im\ > 0, A := fc^, has 
exactly one solution which belongs to L^(R+) (the limit-point case at infinity). 

In this case there is exactly one spectral function p(A) of the selfadjoint operator 
Lq. Denote 

I sin(VAa;) 

(5.1) Lpo{x,X):= -j= — . 

Let h{x) e Lg(IR+), where Lo(I^+) denotes the set of -L^(IR+) functions vanishing 
outside a compact interval (this interval depends on h{x)). Denote 

/■OO 

(5.2) H{X):= h{x)<fio{x,X)dx. 

Jo 

Assume that for every h G Lo(M+) one has: 

/OO 
H^{\) dp{X) = ^ h{x) = 0. 
-OO 

Denote by V the set of nondecreasing functions p{X), of bounded variation, such 
that if pi,p2 GV, u := p\ — p2, and 

(5.4) n:={H{X):hGC^{R+)}, 
where H{X) is given by (5.2), then 

(5.5) I J H'^{X) dv{X) =0 V/i e w| ^ v{X) = 0. 

Theorem 5.1. Spectral functions of the operators Lq, in the limit- point at infinity 
case, belong to V . 

Proof. Let 6 > be arbitrary, / e L^(0, 6), / = if a; > 6. Suppose 

/OO 
H'^{X)dv{X) =Q VhGH. 
-OO 

Denote by I+V and I+W the transformation operators corresponding to potentials 
qi and q2 which generate spetral functions pi and p2, v = Pi — P2- Then 

(5.7) V'o = (/ + V)^i = (7 + W)ip2, 

where V and W are Volterra-type operators. Condition (5.6) implies: 

(5.8) ||(7 + y*)/|| = ||(7 + H^*)/|| V/eL^(0,6), 
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where V* is the adjoint operator and the norm in (5.8) is i^(0, &)-norm. Note that 

(5.9) Vf := r V{x,y)f{y)dy, 

Jo 

and 

(5.10) Vy= fv{y,s)f{y)dy. 

J s 

From (5.8) it follows that 

(5.11) I + V* ^Uil + W*), 

where U is a, unitary operator in the Hilbert space H = L^(0, fe). 

If U is unitary and V, W are Volterra operators then (5.11) implies V = W. 

This is proved in Lemma 5.1 below. If \^ = T'K then ipi{x, A) = f2{x, therefore 
Qi = Q2 and Pi (A) = P2{^)- Here we have used the assumption about being in 
the limit-point at infinity case: this assumption implies that the spectral function 
is uniquely determined by the potential (in the limit-circle case at infinity there 
are many spectral functions corresponding to the given potential). Thus if qi = q2, 
then pi(A) = P2(A). Theorem 5.1 is proved. □ 

Lemma 5.1. Assume that U is unitary and V,W are Volterra operators in H = 

L2(G,6). Then (5.11) implies V = W. 

Proof. From (5.11) one gets I + V = {I + W)U* and, using U*U = I, one gets 

(5.12) {I + V){I + V*) = {I + W){I + W*). 
Denote 

(5.13) (I + V)-^ =I + Vi, (I + W)-^ =I + Wu 
where Vi, Wi are Volterra operators. From (5.12) one gets: 

(5.14) {I + V*){I + Wl) = {I + Vi){I + W), 
or 

(5.15) V* + + V*W^ = Vi + W + ViW. 

Since the left-hand side in (5.15) is a Voletcrra operator of the type (5.10) while 
the right-hand side is a Volterra operator of the type (5.9), they can be equal only 
if each equals zero: 

(5.16) V* + + V*W^ = 0, 
and 

(5.17) Vi+W + ViW = 0. 
Prom (5.17) one gets 

(5.18) Viil + W) = -W, 

or [(/ + y)-i - /](/ + W) = -W. Thus (/ + V)-\I + W) = IandV = Was 
claimed. Lemma 5.1 is proved. □ 



20 



A.G. RAMM 



The inverse spectral problem consists of finding q{x) given p{X). The uniqueness 
of the solution to this problem was proved by Marchenko iQ while the reconstruc- 
tion algorithm was given by Gelfand and Levitan [|l) (see also jr^). 

Let us prove first the uniqueness theorem of Marchenko following |Q . In this 
theorem there is no need to assume that Lq is in the limit-point at infinity case: 
if it is not, the spectral function determines the potential uniquely also, but the 
potential does not determine the spectral function uniquely. 

Theorem 5.2. The spectral function determines q{x) uniquely. 

Proof. If qi and (72 have the same spectral function p(A) then 

/oo j*00 
|Fi(A)pdp(A) = / \F^{X)\^dp{\)^\\gr 
-00 J —00 

for any / G iv^(0, 6), h < oo, where 

(5.20) F,{X):= [ fix)^,{x,X)dx, j = l,2, 

Jq 

the function ipj(x,X) solves equation (1.1) with q = qj, and — A, satisfies first 
two conditions (1.4), and 

(5.21) g:={I + K*)f, 
where I + K is the transformation operator: 

(5.22) ^2 = {I + K)^i=Vi+ I K{x,y)My,X)dy. 
Note that 

(5.23) F2(A)= / f{x)(I + K)^idx^ / g{x)(pi{x,X) dx. 

Jo Jo 

From (5.19) it follows that 

(5.24) 11/11 = ||(/ + if*)/|| \ff eL'{0,b):=H. 

Since Range (/ + A'*) = H, equation (5.24) implies that I + K is unitary (an 
isometry whose range is the whole space H). Thus 

(5.25) I + K = {I + K*)-^ = I + T*, 

where T* is a Volterra operator of the type (5.10). 

Therefore K = T* and this implies K = T* = 0. Therefore (pi — (p2 and qi = q2- 
Theorem 5.2 is proved. □ 

Let dpj{X),j = 1,2, be the spectral functions corresponding to the operators Lq.. 
Assume that dpi (A) — cdp2{X), where c > is a constant. The above argument 
can be used with a minor change to prove that this assumption implies: c = 1 and 
Qi ~ <l2- Indeed, the above assumption implies unitarity of the operator ^{I + K). 
Therefore c{I + K) = I + T*. Thus c = 1 and AT = T* = 0, as in the proof of 
Theorem 5.2. Here we have used a simple claim: 

If bl + Q ~ 0, where b = const and Q is a linear compact operator in H , then 
6 = andQ = 0. 

To prove this claim, take an arbitrary orthonormal basis {«„} of the Hilbert 
space H. Then ^ as n ^ oo since Q is compact. Note that ||u„|| = 1, 
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SO 5 = = IIQi^nll as n oo. Therefore b — and consequently Q — 0. 

The claim is proved. 

The Gelfand-Levitan (GL) reconstruction procedure is: 

(5.26) p(A) ^ L{x, v) => K{x, y) => qix). 

Here 

/oo 
tpo{x,X)ipQ{y,\)d(T{\), da = dp - dpo, 
-OC 

A>0, 
A < 0. 



(5.28) dpQ 



Compare (5.28) and (2.5) and conclude that po is the spectral function corre- 

ith q{x) - 
1.10)) 

sin(a;VA) 



sponding to the Dirichlet operator £q — — ^ + q{x) in L'^(R^) with q{x) = 
The function K{x,y) defines the transformation operator (cf. (1.10)) 



(5.29) ip{x,X) ^ ipo{x,X) + J K{x,y)ipQ{y,X)dy, ipo . ^ ' 

where ip solves (1.1) with fc^ = A and satisifics first two conditions (1.4). 

One can prove (see [|l^ , ) , that K and L are related by the Gelfand-Levitan 
equation: 

(5.30) K{x,y) + L{x,y)+ [ K{x,t)L{t,y) dt ^ 0, 0<y<x. 

Jo 

Let us assume that the data p(A) generate the kernel L{x, y) (by formula (5.27)) 
such that equation (5.30) is a Fredholm-type equation in L^{0,x) for K{x,y) for 
any fixed x > 0. 

Then, one can prove that assumption (5.3) implies the unique solvability of the 
equation (5.30) for K{x,y) in the space L'^{0,x). 
Indeed, the homogeneous equation (5.30) 

(5.31) h{v)+ ( L{t,y)h{t)dt^O, 0<y<x, 

Jo 

implies h — it (5.3) holds. To see this, multiply (5.21) by h and integrate over 
{0,x) (assuming without loss of generaity that h = h, since the kernel L{t,y) is 
real- valued) . The result is: 

0=\\hr+ r \H{X)\^dp{X)- r \H{X)\^dp,{X), 

J —oo J —oo 

or, by Parseval's equality, 

poo 

(5.32) / \H{X)fdp{X) = 0. 



From (5.32) and (5.3) it follows that h{y) — 0. Therefore, by Fredholm's alternative, 
equation (5.30) is uniquely solvable. 
If K{x, y) is its solution, then 

, , dKix.x) 
(5.33) q{x) = 2 

If K{x,x) is a (7™+^ function then q{x) is a C™-function. 
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One has to prove that potential (5.22) generate the spectral function p{X) with 
which we started the inversion procedure (5.26). 
We want to prove more, namely the following: 

Theorem 5.3. Each step in the diagram (5.26) is invertible, so 

(5.34) p<^ L<^ K ^q. 

Proof. 1) Step p L is done by formula (5.27). 

Let us prove L ^ p. Assume there are pi and p2 corresponding to the same 
L(x,y). Then 

/oo 
ipo{xTX)ipo{y,\) dv Vx,?/eM. 
-oo 

Therefore 

/•oo 

(5.36) = / H^{\)dv V/ieC^(M+). 



By Theorem 5.1 relation (5.36) implies v{X) = 0, so pi ~ p2. 

2) Step L ^ K is done by solving equation (5.30) for K{x,y). The unique 
solvability of this equation for K{x,y) has been proved below formula (5.30). 

Let up prove K ^ L. From (5.27) one gets 

Tf ^ L{x + y) - L{x ~ y) /■°° 1 - cos(x\/A) 

(5.37) L{x,y) = , L{x) -.^ J da{X). 

Let y — X in (5.30) and write (5.30) as 

(5.38) L{2x)+ [ K{x,t)[L{x + t)-L{t-x)]dt^-2K{x,x), x>0. 



Note that L{—x) = L{x). Thus (5.38) can be written as: 
(5.39) 

L{2x)+ K{x,s - x)L{s)ds - K{x, x - s)L{s)ds = -2K{x, x), x > 0. 

Jx Jo 

This is a Volterra integral equation for L{s). Since it is uniquely solvable, L[s) 
is uniquely recovered from K(x, y) and the step K ^ L is done. 
3) Step if => g is done by equation (5.33). 

The converse step q ^ K is done by solving the Goursat problem: 

(5.40) K,,~q{x)K ^Kyy 0<y<x, 

(5.41) K{x,x) = - / q{t)dt, K{x,0) = 0. 

2 Jo 

One can prove that any twice differentiable solution to (5.30) solves (5.40)-(5.41) 
with q{x) given by (5.33). The Goursat problem (5.40)-(5.41) is known to have a 
unique solution. Problem (5.40)-(5.41) is equivalent to a Volterra equation (|]l^, 
0)- 

Namely if ^ = x + y, rj ^ x - y, K{x,y) := B{(_,ri), then (5.40)-(5.41) take the 
form 

(5.42) B^,-, = \qi^)B{tv), B{^,0) = ^ J%{t) dt, i?(e,O=0. 
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Therefore 

(5.43) B{^, ^) = lj\(^^jdt+\j^^ dt q {^-±^^ Bit, s) ds. 
This Volterra equation is uniquely solvable for B{^,ri). 

Theorem 5.3 is proved. □ 

5.2. Inverse scattering problem on the half-line. This problem consists of 
finding q{x) given the data (2.10). Theorem 2.2 guarantees the uniqueness of the 
solution of this inverse problem in the class Li.i :— Li.i(M+) of the potentials. 

The characterization of the scattering data (2.10) is known (|Q, ||l^), that is, 
necessary and sufficient conditions on S for S to be the scattering data correspond- 
ing to a q{x) G il l. We state the result without proof. A proof can be found in 
|T7[ . A different but equivalent version of the result is given in [ p^ . 

Theorem 5.4. For the data (2.10) to he the scattering data corresponding to a 
q G ii.i it is necessary and sufficient that the following conditions hold: 

(5.44) i) ind S{k) = -K <0, k^2J or /t = 2 J + 1, 

(5.45) a) kj > 0, sj > 0, l<j<J, 

(5.46) iii) 'S{k) = S{-k) = S-\k), S{oo) = l, keR, 

(5-47) iv) \\F{x)\\l^^u^) + ||F(a;)||ii(R^) + \\xF' {x)\\l^j^^^ < oo. 

Here k = 2./ + 1 if /(O) ^ and k ^ 2J if /(O) ^ 0, and 

(5.48) Fix) := — / [1 - S{k)y''^dk + V s^e"'^'^". 

271- J-oo 



The following estimates are useful (see |J, p.209, [O, |36 , p.569 ) 



\F{2x) + A{x,x)\ < c \qix)\dx, \F{2x)\ < c \q{x)\dx, 

J X J X 

\F'(2x)~'^\<c{£ \q{x)\dx)\ 

where c > is a constant. The Marchenko inversion procedure for finding q{x) 
from iS is described by the following diagram 

(5.49) S => F{x) => A{x, y) => q{x). 

The step 5 =^ is done by formula (5.48). 

The step F ^ A is done by solving the Marchenko equation for A{x, y): 



(5.50) A{x,y)+F{x + y)+ A{x,t)F{t + y) dt = 0, y > x > 0. 

J X 

The step A ^ q is done by the formula 

(5.51) 9(2;) = --'^^^'''''^ 



dx 

It is important to check that the potential q{x) obtained by the scheme (5.49) 
generates the same data S with which we started the inversion scheme (5.49). 
Assuming q G Li^i we prove: 
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Theorem 5.5. Each step of the diagram (5.4-9) is invertible: 

(5.52) S ^ F ^ A<=^ q. 

Proof. 1. The step S ^ F is done by formula (5.48) as we have already mentioned. 

The step F =^ S is done by finding kj, Sj and J from the asymptotics of the 
function (5.48) as a; — > — oo. As a result, one finds the function 

,7 

(5.53) Fd(x) :=^Sje-'^'^^ 

i=i 

If F{x) and Fd{x) are known, then the function 

1 r°° 

(5.54) Fs{x):=— [1 - S {k)^'''^ dk 

is known. Now the fmiction S{k) can be found by the formula 

/oo 
Fs[x)e-'^''dx. 
-oo 

So the step F 5 is done. 

2. The step F A is done by solving equation (5.50) for A{x^y). This step is 
discussed in the literature in detail, (see 0, 0). If q G Li^i (actually a weaker 
condition ^[^(x)!^^ < cx) is used in the half-line scattering theory), then one 
proves that conditions i) - iv) of Theorem 5.4 are satisified, that the operator 

/•oo 

(5.56) r/:= / F{y + t)f{t)dt, y>x>0, 

J X 

is compact in L^{x, oo) and in L^{x, oo) for any fixed a; > 0, and the homogeneous 
version of equation (5.50): 

(5.57) f + Tf = 0, y>x>0, 

has only the trivial solution / = for every a; > 0. Thus, by the Fredholm 
alternative, equation (5.50) is uniquely solvable in L'^(x,oo) and in L^(x,oo). The 
step F ^ A is done. 

Consider the step A{x,y) ^ F{x). Define 

'^(0,2;), 7/>0, 
0, y<0. 

The function A(y) determines uniquely f{k) by the formula: 



(5.58) Aiy) 



(5.59) f{k) = 1 + / A{y)e'''ydy 



and consequently it determines the numbers ikj as the only zeros of f{k) in €+, 
the number J of these zeros, and S{k) = To find F{x), one has to find 

Sj. Formula (2.12) allows one to calculate Sj if f{k) and f'{0,ikj) are known. 
To find f'{Q,ikj), use formula (2.26) and put k — ikj in (2.26). Since A{x,y) 
is known tor y > x > 0, formula (2.26) allows one to calculate f'{0,ikj). Thus 
S{k),kj,Sj,l < j < J, are found and F{x) can be calculated by formula (5.48). 
Step A^ F is done. 

The above argument proves that the knowledge of two functions A{0, y) and 
Ax{Q,y) for all 7/ > determines q{x) G Li_i(R_|_) uniquely. 
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Note that: 

a) we have used the foUowing scheme A ^ S ^ F in order to get the imphcation 
A^ F, 

and 

b) since {F{x),x > 0} ^ A and A {F{x), —oo < x < oo}, we have proved 
also the following non-trivial implication {F{x),x > 0} =J> oo < a: < oo}. 

3. Step A ^ q is done by formula (5.51). The converse step q ^ A is done by 
solving the Goursat problem: 

(5.60) A^^~q{x)A^ Ayy, y>x>0, 



1 f°° 

(5.61) A{x,x) = - / q{t)dt, A(a;, y) — > as x + y ^ oo. 

Problem (5.60)-(5.61) is equivalent to a Volterra integral equation for A{x,y) (see 
p.253]). 



1 



2 



OO 



(5.62) A{x,y) = - q{t)dt+ ds dtq{s - t)A{s - t, s + t). 



One can prove that any twice differentiable solution to (5.50) solves (5.60)-(5.61) 
with q{x) given by (5.51). 

A proof can be found in Q, |n| and |l|]. 

Theorem 5.5 is proved. □ 



Remark 5.1. It follows from Theorem 5.5 that the potential obtained by the scheme 
(5.49) generates the scattering data S with which the inversion procedure (5.49) 
started. 

Similarly, Theorem 5.3 shows that the potential obtained by the scheme (5.26) 
generates the spectral function p(A) with which the inversion procedure (5.26) 
started. 

The last conclusion one can obtain only because of the assumption that q{x) is 
such that the limit-point case at infinity is valid. 

If this is not the case then there are many spectral function corresponding to a 
given q{x), so one cannot claim that the p(A) with which we started is the (unique) 
spectral function which is generated by q{x), it is just one of many such spectral 
functions. 

Remark 5.2. In [^of the following new equation is derived: 

/•oo 

(5.63) F{y)+A{y)+ A{t)F{t + y)dt = A{-y), -00 < y < 00, 

Jo 

which generalizes the usual equation (5.50) at x — 0: 

1*00 

(5.64) F{y)+A{y)+ A{t)F{t + y)dt ^ 0, y > 0. 

Jo 

Since A{—y) = for y > (see (5.58)), equation (5.64) follows from (5.63) for 
y > 0. For y = equation (5.64) follows from (5.63) by taking y — > +0 and using 
(5.58). 
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Let US prove that equations (5.63) and (5.64) o-'"'^ equivalent. Note that equation 

(5.64) is uniquely solvable if the data (2.10) correspond to a g € i^i,!- Since any 
solution of (5.63) in L^(]R+) solves (5.64), any solution to (5.63) equals to the unique 
solution A{y) of (5.64) for y > 0. 

Since we are looking for the solution A{y) of (5.63) such that A = for y < 
(see (5.58)) one needs only to check that (5.63) is satisfied by the unique solution 
of (5.64). 

Lemma 5.2. Equation (5.63) and (5.64) '^''c equivalent in L^{Rj^). 

Proof. Clearly, every i^(R+) solution to (5.63) solves (5.64). Let us prove the 
converse. Let A{y) e iy^(M_|_) solve (5.64). Define 

/•OO 

(5.65) /(fc) := 1 + / A{y)e'^ydy := 1 + A{k). 



We wish to prove that A{y) solves equation (5.63). Take the Fourier transform 
of (5.63) in the sense of distributions. From (5.48) one gets 

/OO J 
F{x)e'^''dx = 1 - 5(-e) + + ikj), 

3 = 1 

and from (5.63) one obtains: 

(5.67) m + m + M-om = i(-e)- 

Add 1 to both sides of (5.67) and use (5.65) to get 

(5.68) fio+mfi-o = fho- 

From (5.66) and (5.68) one gets: 
(5.69) 

J 1 J 

fiO = S{-0 - 27r ^ ,s,<5(e + ik,) /(-O = /(O - 27r/(-e) ^ s,S{C + ik,) = /(O. 

j=l J 3 = 1 

Equation (5.69) is equivalent to (5.63) since all the transformations which led from 
(5.63) to (5.69) are invertible. Thus, equations (5.63) and (5.69) hold (or fail to 
hold) simultaneously. Equation (5.69) clearly holds because 

3 J 

(5.70) /(-c) 'Ai + - E ^^f^'^i) = 0' 

3=1 3=1 

since ikj are zeros of /(fc). 

Lemma 5.2 is proved. □ 

The results and proofs in this section are partly new and partly are based on the 
results in and 0. 

6. Inverse problem for an inhomogeneous Schrodinger equation 
Consider the problem 

(6.1) u" + k^u — q[x)u = —6{x), — OO < x < oo, 

du 

(6.2) ——7 — iku — > as \x\ — > oo. 

d\x\ 
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Assume 

(6.3) q^q, g = for |a;| > 1, geL^[-l,l]. 
Suppose the data 

(6.4) {u{-l,k),u{l,k)}Mk>o 

are given. 

The inverse problem (IP6) is: 
Given data (6.4), findq{x). 
Let us also assume that 

(A): The operator Lq = —-^ + q{x) in L?iW) has no negative eigenvalues. 

This is so if, for example, q{x) > 0. 

The results of this section are taken from |^ 

Theorem 6.1. // (6.3) and (A) hold then the data (6.4) determine q{x) uniquely. 

Proof. The solution to (6.1)-(6.2) is 



(6.5) 



''Mfix^k), x>0, 



(tti 9ix,k), x<0, 



lf.9] 

m 

If.g] 

where / = /+(x,/c), g = g-{x,k), g{k) := g^iO,k), f{k) := f{0,k), [f,g] := 
fg' — f'g = ~2ika{k), a{k) is defined in (2.53), / is defined in (1.3) and g is defined 
in (1.5). 

The functions 

«<-'^^)^ W 

are the data (6.4). 

Since q — Q when x ^ [—1, 1], condition (6.2) implies /(I, k) = e''^, so one knows 

(6.7) '^i^^^^^fH' ^^(^)^={fy' ^^>o- 

From (6.7), (2.49) and (2.50) one derives 
(6.8) 

a{k)hi{k) = -b{-k)f{k) + a{k)f{-k) = -h{-k)h2(k)a{k) + h2{-k)a{~-k)a{k), 
and 

(6.9) a(k)h2{k) = b{k)a(k)hi{k) + a{k)hi{-k)a{-k). 
From (6.8) and (6.9) one gets 

(6.10) -b{-k)h2{k) + ai-k)h2i-k) = hi{k), 
and 

(6.11) b{k)hi{k) + a{~k)hi{-k) = h2{k). 
Eliminate b{—k) from (6.10) and (6.11) to get 

(6.12) a{k)=m{k)a{-k)+n{k), Vfc G R, 
where 
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Problem (6.12) is a Riemann problem for the pair {a{k),a{—k)}, the function 

a{k) is analytic in C+ := {fc : fc € C, Imk > 0} and a(— fc) is analytic in C_. The 
functions o(fc) and a{—k) tend to one as k tends to infinity in C+ and, respectively, 
in C_, see equation (2.55). 

The function a{k) has finitely many simple zeros at the points ikj, 1 < j < J, 
kj > 0, where —kj are the negative eigenvalues of the operator £ defined by the 
differential expression iu = —u" + q{x)u in L^(R). 

The zeros ikj are the only zeros of a{k) in the upper half-plane k. 

Define 



1 f°° 

ind a{k) := 7— / dlna(A;). 
27ri J_^ 



One has 

(6.14) ind a = J, 

where J is the number of negative eigenvalues of the operator £, and, using (6.14) 

and (6.15), one gets 

(6.15) ind m{k) -2[ind hi{k) + ind h2{k)] = -2[ind g{k) + ind f{k) - 2J]. 

Since t has no negative eigenvalues by the assumption (A), it follows that J = 0. 

In this case ind f{k) = ind g{k) = (see Lemma 1 below), so ind m{k) = 0, and 
a(k) is uniquely recovered from the dat,a as the solution of (6.12) which tends to one 
at infinity. If a{k) is found, then b{k) is uniquely determined by equation (6.11) 
and so the reflection coefficient r(fc) := is found. The reflection coefficient 
determines a compactly supported q{x) uniquely by Theorem 2.4. 

If q{x) is compactly supported, then the reflection coefficient r{k) := is 
meromorphic. Therefore, its values for all > determine uniquely r(k) in the 
whole complex fc-plane as a meromorphic function. The poles of this function in 
the upper half-plane are the numbers ikj = 1, 2, . . . , J. They determine imiqucly 
the numbers kj, I < j < J, which are a part of the standard scattering data {r{k), 
kj, Sj, 1 < j < J}, where Sj are the norming constants. 

Note that if a{ikj) = then b(ikj) ^ 0, otherwise equation (2.49) would imply 
f{x,ikj) = in contradiction to (1.3). 

If r{k) is meromorphic, then the norming constants can be calculated by the for- 
mula Sj = — ^ ^(jfc"') = — «R-eSfc=ifc^. r(fc), where the dot denotes differentiation with 
respect to k, and Res denotes the residue. So, for compactly supported potential 
the values of r{k) for all A; > determine uniquely the standard scatering data, 
that is, the reflection coefficient, the bound states — and the norming constants 
Sj,^ < j < J- These data determine the potential uniquely. 

Theorem 6.1 is proved. □ □ 



Lemma 6.1. If J = then ind f = ind g = 0. 

Proof. We prove ind f = 0. The proof of the equation ind g = is similar. Since 
ind f{k) equals to the number of zeros of f{k) in C+, we have to prove that f{k) 
does not vanish in C+. If f{z) = 0, 2: e C+, then z = ik, k > Q, and — fc^ is an 
eigenvalue of the operator I in 1/^(0, 00) with the boundary condition u(0) = 0. 

From the variational principle one can find the negative eigenvalues of the oper- 
ator i in L^(R+) with the Dirichlet condition at a; = as consequitive minima of 
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the quadratic functional. The minimal eigenvalue is: 

(6.16) inf / [u'^ + q{xy]dx:=Ko, ueH\R+), ||u|U2(r^) = 1, 

o 

where ^(M-|-) is the Sobolev space of i?^(M+)-functions satisfying the condition 

u(0) = 0. 

On the other hand, if J = 0, then 



/oo 
lu'"^ + q(x)u^]dx := Ki, ueH^(M.), \\u\\l2( 
-oo 



= 1. 



Since any element u of H ^(K+) can be considered as an element of H^{M.) if 
one extends u to the whole axis by setting u — for a; < 0, it follows from the 
variational definitions (6.16) and (6.17) that ki < kq. Therefore, if J = 0, then 
Ki > and therefore kq > 0. This means that the operator £ on i^(M_|_) with the 
Dirichlet condition at x = has no negative eigenvalues. Therefore f{k) does not 
have zeros in C+, if J = 0. Thus J = implies ind f{k) = 0. 

Lemma 6.1 is proved. □ 

The above argument shows that in general 

(6.18) ind f < J and ind g < J, 
so that (6.15) implies 

(6.19) ind m{k) > 0. 

Therefore the Riemann problem (2.17) is always solvable. It is of interest to 
study the case when assumption (A) does not hold. 

7. Inverse scattering problem with fixed energy data 

7.1. Three-dimensional inverse scattering problem. Property C. The scat- 
tering problem in M"^ consists of finding the scattering solution u := u{x, a, k) from 
the equation 

(7.1) [V^ + A;^ - g(a;)] V = in 
and the radiation condition at infinity: 

(7.2) V = V'o + t^, Vo:=e^'=«•^ a e 5^, 



(7.3) lim / 



dv 
o\x\ 



ds = 0. 



I\s\=r 

Here A; > is fixed, S'^ is the unit sphere, a e 5^ is given. One can write 

(7.4) v = A{a',a,k) as r = |a;| — > oo, — = a! . 

r \r J r 

The coefficient A{a' , a, k) is called the scattering amplitude. In principle, it can 

be measured. We consider its values for a' , a E S'^ and a fixed fc > as the 
scattering data. Below we take fc = 1 without loss of generality. 
Assume that 

(7.5) q&Qa-={(l-q = q, 9 = for |a;| > a, q& LP{Ba)} , 
where a > is an arbitrary large fixed number, Ba = {x : \x\ <a},p> |. 
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It is known (even for much larger class of the potentials q) that problem (7.1)- 
(7.3) has the unique solution. 

Therefore the map q — > A{a' , a) :— Aq(a' , a) is well defined, 

A{a',a) := A{a' , a,k)\k=i. 

(IP7) The inverse scattering problem with fixed-energy data consists of finding 
q{x) G Qa from the scattering data A{a',a) Va', a G . 

Uniqueness of the solution to IP7 for q G Qa (with p = 2) was first announced in 
|27| and proved in p8| by the method, based on property C for pairs of differential 
operators. The essence of this method is briefly explained below. This method was 
introduced by the author and applied to many inverse problems |27|-||3^, p^ . 

In [|17| and a characterization of the fixed-energy scattering amplitudes is 
given. 

Let {Li,L2} be two linear formal differential expressions, 
(7.6) LjU= V a,„j(a:)a'"u(a:), x G M", n > 1, j = 1,2, 



|m|=0 



where 



Q\m\ 

^" ■= ^-^fn fJJT' +m„. 

OX-^ . . . Xji 

Let 

(7.7) Nj := Nj{D) := {w : LjW = in £> C i?"} 

where D is an abitrary fixed bounded domain and the equation in (7.7) is under- 
stood in the sense of distributions. 
Suppose that 



(7.8) f{x)wi{x)w2{x)dx = Q, w^eNj, f e L'{D), 

J D 

where Wj G A'^ run through such subsets of Nj, j ~ 1, 2, that the products wiW2 G 
L^{D), and / G L^{D) is an arbitrary fixed function. 

Definition 7.1. The pair {Li,L2} has property C if (7.8) implies f{x) = 0, that 
is, the set {wi,W2}\iwjeNj, wiWiEL^iD) is complete in L'^{D). 

In | |3l| a necessary and sufficient condition is found for a pair {Li,L2} with 
constant coefficients, Omjix) = Omj = const, to have property C (see also [pT[). 

In 1^ it is proved that the pair {Li,L2} with Lj = — + q{x), qj G Qa, has 
property C. 

The basic idea of the proof of the uniqueness theorem for inverse scattering 



problem with fixed-energy data, introduced in |27[ , presented in detail in |28 , and 
developed in [^-|3^, is simple. Assume that there are two potentials, qi 

and q2 in Qa which generate the same scattering data, that is, Ai — A2, where 
A, :=^,^.(a',a), j = l,2. 
We prove that 0, p. 67] 



(7.9) ~A'k{Ai-A2)^ / [qi{x)—q2{x)]'4'i{x,a)4>2{x,—a')dx ya,a'eS'^, 

JBa 

where ijjj{x,a) is the scattering solution corresponding to qj, j — 1, 2. 
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If Ai — A2, then (7.9) yields an orthogonality relation: 

(7.10) / p{x)ipi{x,a)ip2{x,l3)dx ^0 Va,/3eS'^, p(x) := gi - 92- 

JD 

Next we prove p. 45] that 

(7.11) span„gs2{V'(a;,a)} is dense in L^{D) in Nj{D) n H^{D), 
where H"'^{D) is the Sobolev space. Thus (7.10) implies 

(7.12) / p{x)wi{x)w2{x) dx^O \fwj e Nj{D) n H'\D). 

JD 

Finally, by property C for a pair {ii, L2}, Lj = — V^+qj, qj G Qa, one concludes 
from (7.11) that p{x) — 0, i.e. qi — q2. We have obtained 

Theorem 7.1 (Ramm [|8|, Q). The data A{a',a) Va', a £ determine q e 
Qa uniquely. 

This is the uniqueness theorem for the solution to (IP7). In fact this theorem is 
proved in [T^ in a stonger form: the data Ai{a',a) and A2{a',a) are asumed to 
be equal not for all a' ,a £ but only on a set x S2, where 5| is an arbitrary 
small open subset of S*^. 

In the stability estimates for the solution to (IP7) with noisy data are ob- 
tained and an algorithm for finding such a solution is proposed. 

The noisy data is an arbitrary function Ag(a\a), not necessarily a scattering 
amplitude, such that 

(7.13) sup \A{a',a) - A^{a',a)\ < e. 

Given Ae(a', a), an algorithm for computing a quantity q^ is proposed in [ p4[ , such 
that 

(7.14) supl^.-.lOKcMl!^. 

where c > is a constant depending on the potential but not on e, 

(7.15) m / e^^--qix)dx. 



The constant c in (7.14) can be chosen uniformly for all potentials q E Qa which 
belong to a compact set in L^{Ba)- 

The right-hand side of (7.14) tends to zero as £ 0, but very slowly. 

The author thinks that the rate (7.14) cannot be improved for the class Qa, but 
this is not proved. 

However, in [Q| an example of two spherically symmetric piecewise-constant po- 
tentials q{r) is constructed such that |gi — is of order 1, maximal value of each 
of the potentials q is of order 1, the two potentials are quite different but they 

generate the set of the fixed-energy ik = 1) phase shifts {5^^ \ , j ~ 1,2, 

I J ^=0,1, 2,... 

such that 



(7.16) (5f^=(5f\ 0<^<4, 



< lO-\ V£ > 5. 



In this example e ^ 10 ^, (In | lne|)^ ~ 2.59, | Inej ~ 5, so the right-hand side of 
(7.14) is of order 1 if one assumes c to be of order 1. 
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Our point is: there are examples in which the left and the right sides of estimate 
(7.14) are of the same order of magnitude. Therefore estimate (7.14) is sharp. 

7.2. Approximate inversion of fixed-energy phase shifts. Let us recall that 
q{x) ~ q{r), r = \x\ if and only if 

(7.17) A{a',a,k) ^ A{a' ■ a,k). 

It was well known for a long time that if q = q{r) then (7.17) holds. The converse 
was proved relatively recently in Q (see also [p7[). 
li q = q{r) then 

00 

(7.18) A(a', a)^J2 ^eYi{a')Yiiaj, 

where Yi are orthonormal in L^{S^) spherical harmonics, Yi = Y^m, —£< m < £, 
summation with respect to m is understood in (7.18) but not written for brevity, 
the numbers Ai are related to the phase shifts Si by the formula 

(7.19) Ae = ine'^' sm(Si) (k = 1), 
and the S'- matrix is related to A{a' , a, k) by the formula 

S = I-^A. 

If g = q{r), r — \x\, then the scattering solution ^-^{x^ a, k) can be written as 



(7.20) V(^,a,A:)= V^/^M!l^i.(^0)^^ ^0 ^ 

^ — ' k r 



e=o 



The function ipi{r, k) solves (uniquely) the equation 

/>oo 

(7.21) ^pi{r, k) ui{kr) - gt{r, p)q{p)-)pe{p, k) dp, 

Jo 

where 



(7.22) ui{kr) := ^ _ J,^,/2(fcr), vi := ^ —Ni+^/^ikr). 

Here Je and Ni are the Bessel and Neumann functions, and 

{ Voi{kp)foe{kr) ^ > „ 
Fot(k) ' ^ ^ ^' 

(7.24) Foiik) = 
ipoi and foe solve the equation 

(7.25) + ^Vo£ - ^^^t^V o^ - 0, 
and are defined by the conditions 

(7.26) /o£ ~ e'''^ as r ^ +00, 
so 

(7.27) foeikr) =ie'^{u(; + ive), 
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where ui and vi are defined in (7.22), and 

ui{kr) 



(7.27') ipM{kr) 



In p9[ | an approximate method was proposed recently for finding q{r) given 
{(5f}£=o,i,2,...- In numerical results based on this method arc described. 

In physics one often assumes q{r) known for r > a and then the data {i5£}£=o,i,2,... 
allow one to calculate the data ipe{a),£ = 0, 1, 2, ... , by solving the equation 



(7.28) V" + V-f - ^ ^2 - = 0, r>a 
together with the condition 

(7.29) ipe^ e'^'am{r-^+6e), r ^ +oo, k = 1, 

and assuming q{r) known for r > a. 

Problem (7.28) and (7.29) is a Cauchy problem with Cauchy data at infinity. 
Asymptotic formula (7.29) can be differentiated. If the data ipeia) are calculated 
and ipi{r) for r > a is found, then one uses the equation 



(7.30) i;,{r) = ijf\r) ~ g,{r, p)q{p)Mp) dp, < r < a, 

Jo 

where ge is given in (7.23), 

/•oo 

(7.31) tpf\r):^ue{r)~ gi{r, p)q{p)^p(,{p) dp, 

J a 

and ui(r) is given in (7.22). 
Put r = a in (7.30) and get 

(7.32) / g,{a,p)Mp)<lip)dp = 4''\'^)-Ma)-=bi, 1 = 0,1,2, 



The numbers bi are known. If q{p) is small or £ is large then the following 
approximation is justified: 

(7.33) Mp)'^4°\p)- 
Therefore, an approximation to equation (7.32) is: 

(7.34) / Mp)q{p)dp = b,, £^0,1,2,..., 
Jo 

where 

(7.35) Mp):=giia,p)4"\p). 

The system of functions {fi{p)} is linearly independent. 

Equations (7.34) constitute a moment problem which can be solved numerically 
for qip) (0 p.209], [|o|). 
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8. A UNIQUENESS THEOREM FOR INVERSION OF FIXED-ENERGY PHASE SHIFTS. 

From Theorem 7.1 it follows that ii q = q{r) e Qa then the data 

(8.1) {6i} V^ = 0,l,2,... fc = l 

determine q{r) uniquely |p8|| . 

Suppose a part of the phase shifts is known (this is the case in practice) . 
What part of the phase shifts is sufficient for the unique recovery of q{r) ? 
In this section we answer this question following pi]. Define 



(8.2) C 



to be any subset of nonnegative integers such that condition (8.2) is satisfied. 
For instance, C = {2£}f=o,i,2,... or £ = {2i+ l}f=o,i,2,... will be admissible. 
Our main result is 

Theorem 8.1 (Q). If q{r) G Qa and L satisfies (8.2) then the set of fixed- energy 
phase shifts {6i\yi£c determine q{r) uniquely. 

Let us outline basic steps of the proof. 

Step 1. Derivation of the orthogonality relation: 

If <Zi, 92 G Qa generate the same data then p{r) := qi — q2 satisfies the relation 



(8.3) / p{r)tl:f\r)i)f\r)dr^Q G £. 



Here 'ip'f\r) are defined in (7.20) and correspond to q ~ qj, j = 1,2. Note that 
Ramm's Theorem 7.1 yields the following conclusion: if 

/•a 

(8.4) / p(r)v4^^ (r)V'f ^ (r) dr ^ G £, 



then p{r) = 0, and qi — q2. 

Step 2. Since ijj'f'^ = c'pLp"p(r), where are some constants, relation (8.3) is 
equivalent to 

/•a 

(8.5) / p{r)ifP {r)ipf\r) dr = Q G £. 

Here fPir) is the solution to (7.28) with q — qj, which satisfies the conditions: 

(8.6) ^^^(2ZTT)!!+°^"''"'^' '^^^ 
and 

(8.7) ^« ^|^^«|sin(r-y +5f ) + o(l), r -> +(X3, 

where {Fp^'l 7^ is a certain constant, and 6^ is the fixed-energy (fc = 1) phase 
shift, which does not depend on j by our assumption: = S^p \/l G C. 

We want to derive from (8.5) the relation (7.10). Since q{x) — q{r), r — \x\ in 
this section, relation (7.10) is equivalent to relation (8.12) (see below). The rest of 
this section contains this derivation of (8.12). 
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We prove existence of the transformation kernel K{r, p), independent of which 
sends functions u^{r), defined in (7.22), into (pe{r)'. 

(8.8) Mr)=Uiir)+ f Kir,p)u,ip)%, i^(r,0) = 0. 

Jo r 

Let 

(8.9) 7<':=/|rQ)2^+ir(^+l), 
where T{z) is the gamma function, 

(8.10) H{e) ■.= ^j f p{r)uj{r)dr, 

Jo 

and 

(8.11) hii) rpir)4'\r)^f\r)dr. 

Jo 

We prove that: 

(8.12) {h{£) = Q Vi e C} =^ {h{e) = V£ = 0,l,2,...}. 

If h{£) = 0, = 0, 1, 2, . . . , then (8.4) holds, and, by theorem 7.1, p{r) = 0. 
Thus, Theorem 8.1 follows. 
The main claims to prove are: 

1) Existence of the representation (8.8) and the estimate 



.13) r \K{r,p)\^ <c{r) <oo 

Jo P 



Vr > 0. 



2) Implication(8.12). 

Representation (8.8) was used in the physical hterature (||^, iQ) but, to our 
knowledge, without any proof. Existence of transformation operators with kernels 
depending on £ was proved in the literature p3[ . For our purposes it is important 
to have K{r, p) independent of £. 

Implication (8.12) will be estblished if one checks that h{t) is a holomorphic 
function of £ in the half-plane :—{£:£ ^ Re£ > 0} which belongs to 
N-Class (Nevanlinna class). 

Definition 8.1. A function h{£) holomorphic in belongs to N-class iff 

■ 1 _ re"^ ' 



(8.14) sup / In^ 

0<r<l J-TT 

Here 



h 



1 + re'-'P 



d(f < oo. 



I 1 7i if w > 0, 
:= < 

|0 if w < 0. 

The basic result we need in order to prove (8.12) is the following uniqueness theo- 
rem: 

Proposition 8.1. If h{£) belongs to N-class then (8.12) holds. 
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Proof. This is an immediate consequence of the following: 

Theorem (]|4l|, p.334]/' If h{z) is holomorphic in Di :— {z : |z| < 1, z G C}, 
h(z) is of N- class in Di, that is: 



(8.15) 

and 
(8.16) 

where 

(8.17) 

then h{z) = 0. 



sup / ln+ \h{re''^)\dif) < oo, 

0<r<l J-TT 

/i(z„) = 0, /i = 1,2,3,..., 

OO 

^(1 - |z„|) = oo, 

h=l 



The function z = maps conformally Y[+ onto Di, £ 



1 + 2 



and if h{£) 



ye € C, then /(z) /i(y^) is holomorphic in Di, f{zi) = for ^ e £ and 



Yipjj and 



(8.18) 



El 



1+i 



< 



1 + E 1- 



eec 



ITT 



= 1 



From (8.18) and the above Theorem Proposition 8.1 follows. 

Thus wc need to check that function (8.11) belongs to N-class, that is, 
holds. 

So step 2, will be completed if one proves (8.8), (8.13) and (8.14). 
Assuming (8.8) and (8.13), one proves (8.14) as follows: 

i) First, one checks that (8.14) holds with H{e) in place of h{e). 

ii) Secondly, using (8.8) one writes h{e) as: 



□ 
.14) 



(8.19) 



hit) = H{i)+-fj / [Ki{r,p)+K2{r,p)]u,{p)ut{r^ 



+ 11 



dp 



"'0 



'^'^Ki{r, t)K2{r, s)ui{t)ui{s) . 



Let us now go through i) and ii) in detail. 
In 0, 8.411.8] one finds the formula: 



.20) 



where is defined in (8.9). 

From (8.20) and (8.10) one gets: 



.21) 



H{e)\ < / dr\p{r)\r 



2t+2 



1 2 

{l-t^fe'"-' dt 



< ca^" , e^a + iT, cr > 0. 



One can assume a > 1 without loss of generality. Note that 
.22) ln+ [ah] < ln+ a + ln+ h if a, 6 > 
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Thus (8.21) implies 



(8.23) 



ln+ 



H 



1 - re''^ 



1 + re*' 
< Unci +21na 



dip < 



In^ 



2B.e 



dtp 



Re- 



1 - re"^ 



< line 



21na 



1 + re^'P 
l-r2 



d(p 



_^ 1 + 2rcos(/5 + 
Here we have used the known formula: 

.24) r 2. 



diy9 < I In c| + 47r In a < cx). 



< r < 1. 



1 + 2r cos If + 1 — 

Thus, we have checked that H{£) e A^(n+)> that is (8.14) holds for H{£). 
Consider the first integral, call it in (8.19). One has, using (8.20) and 

(8.13), 

Ihm < r dpp-\\K,{r,p)\ + \K,{r,p)\y+^p' < c{a)a'^ , 
Jo 

(8.25) e = a + iT, a>0. 

Therefore one checks that satisfies (8.14) (with Ji in place oih) as in (8.23). 
The second integral in (8.19), call it hi^), is estimated similarly: one uses (8.20) 
and (8.13) and obtains the following estimate: 

(8.26) 1/2(^)1 < c(a)a2'", i^a + ir, a>0. 

Thus, we have proved that h{i) e iV(n+)- 

To complete the proof one has to derive (8.4) and check (8.13). 

Derivation of (8.4): Subtract from (7.28) with q = qi this equation with q = q2 
and get: 



(8.27) 

where 
(8.28) 



w + w ■ 



-w — qiw = pi/jg 



(2) 



p:=qi-q2, w := ^p'l^\r) - ipf (r) . 

Multiply (8.27) by ■0^^''(»'), integrate over [0, 00), and then by parts on the left, and 
get 



(8.29) (w'ip'^/'> - wV'f ^') 



By the assumption jI^"* = S^^"^ if £ e £, so w and w' vanish at infinity. At r = 
the left-hand side of (8.29) vanishes since 

(8.30) ^//\r) = 0{r'^+^) as r -> 0. 

Thus (8.29) implies (8.4). □ 

Derivation of the representation (8.8) and of the estimate (8.13). One can 
prove ||2^ that the kernel K{r,p) of the transformation operator must solve the 
Goursat-type problem 

(8.31) r'^Krrir, p) - p'Kppir, p) + p - r''q{r) - p']K{r, p) = 0, < p < r. 
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(8.32) i^(r,r) = ^ sq{s) ds := g{r), 

(8.33) K{r,Q) = Q, 

and conversely: the solution to this Goursat-type problem is the kernel of the 
transformation operator (8.8). 

The difficulty in a study of the problem comes from the fact that the coefficients 

in front of the second derivatives degenerate at p = 0, r = 0. 
To overcome this difficulty let us introduce new variables: 

(8.34) ^ = lnr + ln/3, 77 = In r — In p. 
Put 

(8.35) K{r,p):=B{^,n). 
Then (8.31)-(8.33) becomes 

(8.36) B^ri-\Br, + Q{i,ri)B = Q, r? > 0, -00 < ^ < 00, 

(8.37) B{^,Q)=g{ei):=G{0, 

(8.38) B(-oo,r7) = 0, > 0, 
where g{r) is defined in (8.32) and 

(8.39) Q(^, r?) := \ [e«+'' - e^+^q (e^) 
Note that 



(8.40) sup e-5G(0 < c, 

— 00<^<00 

(8.41) sup / \Q{s,ri)\ds<c{A,B), 

0<r,<B J-00 

for any ^ G R and any B > 0, where c{A, i?) > is some constant. 
Let 

(8.42) L(^,7?):=B(^,r?)e-5 
Write (8.36)-(8.38) as 

(8.43) L^^ + Q(^, v)L = 0, > 0, -00 < ^ < 00 

(8.44) L(e,0) = e-iG(O :-6(C); L(-oo,r?)=0, r, > 0. 

Integrate (8.43) with respect to 77 and use (8.44), and then integrate with respect 
to ^ to get: 

(8.45) L = VL + b, VL:=- ds T dtQ{s,t)L{s,t). 

J-00 Jo 

Consider (8.45) in the Banch space X of continous function L{^,r]) defined for 
r? > 0, —00 < ^ < 00, with the norm 

(8.46) ||L|| := \\L\\ab := sup (e-T*|i(s, t)|) < 00, 

0<t<B 
-oo<s<A 
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where 7 = 7(A, B) > is chosen so that the operator F is a contraction mapping 
in X. Let us estimate 



(8.47) 

\\VL\\< sup j ds dt\Q{s,t)\e-^^'>-*^e-'^*\L{s,t)\ 

-oo<{<aJ_oo Jo 



0<?7<B 



< ll^^ll 

0<ri<B 



sup r ds f ' dt (2e'+* + e'+' 
xj<5<aJ_oo Jo ^ 



q[e 2 



-7(77-t) ^ c 

~ 7 



where c > is a constant which depends on A, B, and on r\q{r) \ dr. 

If 7 > c then y is a contraction mapping in X and equation (8.45) has a unique 
solution in X for any — 00 < ^4 < 00 and B > 0. 

Let us now prove that estimate (8.13) holds for the constructed function K{r, p). 

One has 

(8.48) / \K{r,p)\p~^dp = r \L{2\nr - r],T])\e~^ dt] < 00 
Jo Jo 

The last inequality follows from the estimate: 

(8.49) \L{^,V)\ < ce(2+^^)['"'i 
where £1 and S2 > are arbitrarily small numbers, 

(8.50) Mi(0:=y dsp{s), p{s) := - {l + \q {e^)\) 

The proof of Theorem 8.1 is complete when (8.49) is proved. 
Lemma 8.1. Estimate (8.49) holds. 
Proof. From (8.45) one gets 

(8.57) m{^,r]) <c+Wm 

where 

(8.52) CO- sup \b{0\<l: r siq{s)ds, m(C,r?) := |L(C,??)|, 

-oo<{<oo ^ Jo 

and 

(8.53) Wm:= [ ds [ dt p{s + t)m{s,t). 

J -00 Jo 

Without loss of generality we can take co = 1 in (8.51): If (8.49) is derived from 

(8.51) with Co = 1, it will hold for any co > (with a different c in (8.49)). Thus, 
consider (8.51) with co = 1 and solve this inequality by iterations. 

One has 

/? fV fV 

ds p{s + t)dt= p-ii^ + t)dt<r]fj,i{^ + r]). 
-00 Jo Jo 

One can prove by induction that 
(8.55) W-1 < 

hi n\ 
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Therefore (8.57) with cq — 1 impUes 



(8.56) m(^,7?)<l + ^ 

Consider 



ni ni 

h=l 



.57) F{z) := 1 + ^ 



Z 



1^2- 



^ (n!) 

n=l ^ ' 

This is an entire function of order ^ and type 2. 
Thus 

(8.58) |F(z)| <ce(2+^i)l^l*^''. 

From (8.56) and (8.58) estimate (8.49) follows. 

Lemma 8.1 is proved. □ 

Theorem 8.1 is proved. □ 
9. Discussion of the Newton- Sabatier procedure for recovery of qir) 

FROM THE FIXED-ENERGY PHASE SHIFTS 

In 1^ and p6) the following procedure is proposed for inversion of the fixed- 
energy phase shifts for q{r). We take k — \\vl what follows. 

Step 1 . Given {5^ jv^^o i 2 one solves an infinite linear algebraic system ((12.2.7) 
in |) 

00 

(9.1) tan(5£ = M«/(l + tan(5£tan(5^/)a^' 

£'=0 

for constants a^. Here 



(9.2) Mw = 



if |£ — ^'1 is even or zero 



,^._V_<,^._V if l^-^'|isodd. 



Assuming that (9.1) is solvable and ag are found, one calculates (see formula (12.2. 
in I) 

(9.3) 



/ 2 \ I '^'^^ (1 tan^ 6i) I 
ci = ae{l + tan Si) il M + 2 Mu>ai>{ta.nSe> - ta.n.Si) > 

I e'=o J 

Step 2 . If C£ are found, one solves the equation for K{r, p) (see formula (12.1.12) 
ini) 



-1 



(9.4) Kir,s)^fir,3)- K{r,t)f{t,s)t-' dt, 

Jo 

where 

00 

(9.5) /(r, s) := ^ ceui{r)ui{s), 

e=o 

and Ui{r) are defined in (7.22). 
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Note that in this section the notations from are used and by this reason the 
kernel K{r,t) in formulas (9.4) and (9.6) differs by sign from the kernel K{r,p) in 
formula (8.8). This explains the minus sign in formula (9.6). 

Assuming that (9.4) is solvable for all r > 0, one calculates 

/ , , , 2d K(r,r) 

(9.6) q{r) = — - 

r ar r 

Note that if (9.4) is not solvable for some r = r > 0, then the procedure breaks 
down because the potential (9.6) is no longer locally integrable in (0, oo). In jl6j it is 
argued that for sufficiently small a equation (9.4) is uniquely solvable by iterations 
for all < r < a, but no discussion of the global solvability, that is, solvability for 
all r > 0, is given. It is assumed in and that the sequence {ci\ in (9.3) does 
not grow fast. In ([^ (12.2.2)]) the following is assumed: 

oo 

(9.7) ^|Q|r2<oo. 

1=1 

Under this assumption, and also under much weaker assumption 
(9.7') IqI < ce"' 

for some 6 > arbitrary large fixed, one can prove that the kernel (9.5) is an entire 
function of r and s. This follows from the known asymptotics of ui{r) as ^ — > oo: 

(9.8) ..(.) = /f(^) ^ ^[1 + 0(1)], ^^+00. 

Thus, equation (9.4) is a Frcdholm-type equation with kernel which is an entire 
function of r and s. Since K{r, 0) = and /(r, s) = /(s, r) = at s = 0, equation 
(9.4) is a Fredholm equation in the space of continuous functions C(0,r) for any 
r > 0. 

If (9.4) is uniquely solvable for all r > 0, then one can prove the following: 

Claim. K{r, s) is an analytic function of r and s in a neighborhood A of the 
positive semiaxis (0, 00) on the complex plane of the variables r and s. 

This claim is proved below, at the end of this section. 
Therefore the potential (9.6) has to have the following: 

Property P: q{r) is an analytic function in A with a possible simple pole atr = 0. 

Most of the potentials do not have this property. Therefore, if one takes any 
potential which does not have property P, for example, a compactly supported 
potential q{r), and if it will be possible to carry through the Newton-Sabatier 
procedure, that is, (9.1) will be solvable for ae and generate ci by formula (9.3) 
such that (9.7) or (9.7') hold, and (9.4) will be uniquely solvable for all r > Q then 
the potential (9.6), which this procedure yields, cannot coincide with the potential 
with which we started. 

An important open question is: assuming that the Newton-Sabatier procedure 
can be carried through, is it true that the reconstructed potential (9.6) generates 
the scattering data, that is, the set of the fixed-energy phase shifts {5i\ with which 
we started. 

In H, pp. 203-205] it is claimed that this is the case. But the arguments in |^ are 
not convincing. In particular, the author was not able to verify equation (12.3.12) 
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in 1^ and in the argument on p. 205 it is not clear why A'^ and S'^ satisfy the same 
equation (12.2.5) as the orginal Ag and 5^. 

In fact, it is claimed in (12.5.6)] that 5'^^ = 0( j), while it is known |37) that if 
q{r) = for r > a and q{r) does not change sign in some interval (a — £, a), q ^ 
if r G (a — e, a), where e > is arbitrarily small fixed number, then 

/If 

(9.9) lim 

e^oo \ e 

It follows from (9.9) that for the above potentials the phase shifts decay very 
fast as ^ — > cx), much faster that j. Therefore 6'^ decaying at the rate j cannot be 
equal to Si, as it is claimed in [|| p. 205], because formula (9.9) implies: 



ea[l + o(l)] 



21 



21 

It is claimed in p. 105, that the Newton-Sabatier procedure leads to "one 
(only one) potential which decreases faster that 5 " and yields the original phase 
shifts. However, if one starts with a compactly supported integrable potential (or 
any other rapidly decaying potential which does not have property P and belongs 
to the class of the potentials for which the uniqueness of the solution to inverse 
scattering problem with fixed-energy data is established), then the Newton-Sabatier 
procedure will not lead to this potential as is proved in this section. Therefore, 
either the potential which the Newton-Sabatier procedure yields does not produce 
the original phase shifts, or there are at least two potentials which produce the 
same phase shifts. 

A more detailed analysis of the Newton-Sabatier procedure is given by the author 
in|. 

Proof of the claim. The idea of the proof is to consider r in (9.4) as a parameter 
and to reduce (9.4) to a Fredholm-type equation with constant integration limits 
and kernel depending on the parameter r. Let t = rr, s = ra, 

fmn\ K{r,rT) f{rT,r(j) f{r,(j) , . 

(9.10) —b{T;r), ■.= a{a,T,r), —g{a;r) 

T rra a 

Then (9.4) can be written as: 



1 



(9.11) b{a;r) — g{r;a) — / a{a,T;r)b{T;r) dr. 

Jo 

Equation (9.11) is equivalent to (9.4), it is a Fredholm-type equation with kernel 
a{a,T;r) which is an entire function of a and of the parameter r. The free term 
g(r, a) is an entire function of r and a. This equation is uniquely solvable for all 
r > by the assumption. Therefore its solution b{a; r) is an analytic function of r 
in a neighborhood of any point r > 0, and it is an entire function of a Thus 
K{r,r) = &(l,r) is an analytic function of r in a a neighborhood of the positive 
semiaxis (0, oo). 

□ 

10. Reduction of some inverse problems to an overdetermined 

Cauchy problem 

Consider, for example, the classical problem of finding q{x) from the knowledge 
of two spectra. Let u solve (1.1) on the interval [0,1] and satisfy the boundary 



PROPERTY C FOR ODE AND APPLICATIONS TO INVERSE PROBLEMS. 



43 



conditions 

(10.1) m(0) = u(1) = 0, 

and let the corresponding eigenvalues fc^ :— A„, n = 1, 2, ... , be given. If 

(10.3) m(0) = u'(l) + /im(1) = 0, 

then the corresponding eigenvalues are n = 1, 2, . . . . 

The inverse problem (IPIO) is: 

Given the two spectra {A„} U {fin}, n — 1, 2, 3, ... , find q{x) (and h in (10.3)). 
Let us reduce this problem to an overdetermined Cauchy problem. Let 



(10.4) u{x,k) 



sin(fca;) 



+ / K{x,y) 



si'n(ky) 
k 



dy (I + K) 



sin(fca;) 



solve (1.1). Then (10.1) and (10.2) imply 
sin \/A^ 



(10.5) 
and 







K(l,y) 7= — dy, n=l,2, 



(10.6) = cosVM^+-ft^(l,l) 







sm ^/i: 



'An 



+ / K4l,y) 



sm{^/ji:;^y) 



It is known Q that 

(10.7) 

and 



A„ = (nvr)^ + Co + o(l), 



oo, 



dy, n=l,2, 



Mn = Ti"^ ( ?i + 77 ) +ci+o(l), n oo, 



(10.8) 

where co and ci can be calculated explicitly, they are proportional to q{x) dx. 
Therefore, 

(10.9) 

1 



l + O 



A„ = nn 













vr - 






, n - 


-!■ OO 













It is known ||l^ that if (10.9) holds then each of the systems of functions: 

(10.10) Wm{y^x)}n=l,2,..., {sin{y/Ji^x)}n=l,2,... 

is complete in L^[Q, 1]. 

Therefore equation (10.5) determines uniquely {K{l,y)}o<y<i and can be used 
for an efficient numerical procedure for finding K{l,y) given the set {Xn}n=i,2,... 
Note that the system {sin(v'A^a;)}„=i_2.... forms a Riesz basis of L^[0, 1] since 
the operator I + K, defined by (10.4) is boundedly invertible and the system 
{u{x, •\/A^)}n=i.2,... forms an orthornormal basis of L^{0, 1). 

Equation (10.6) determines uniquely {Kx{l, y)}o<y<i if {fJ'n}n=i,2,... sue known. 
Indeed, the argument is the same as above. The constant i4r(l,l) is uniquely 
determined by the data {nn}n=i,2,... because, by formula (5.41), 



1 

if(l,l)-- q(x)dx. 



(10.11) 

We have arrived at the following overdetermined Cauchy problem: 
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Given the Cauchy data 

(10.12) {K{l,y),K^il,y)}o<y<i 

and the equations (5.40) - (5.41), find q{x). 

It is easy to derive [p2l eq. (4.36)] for the unknown vector function 



(10.13) U 



q{x) 
K{x,y) 



the foUowing equation 

(10.14) U = W{U) + h, 

where 

Dxy is the region bounded by the straight hues on the (s, t) plane: s — I, t — y = 
s — X and t — y = —(s — x), and 

(10.16) h := (^^ 

(10.17) f{x) := 2[Ky{l, 2x^1)+ K^l, 2x - 1)], 

(10.18, -ir'"/«M).«. 

^ ^ Jy+x-l 

Note that / and g are computable from data (10.12), and (10.14) is a nonhnear 
equation for q{x) and K{x,y). 
Consider the iterative process: 

(10.19) C/„+i = W{Un) + h, Uo = h. 
Assume that 

(10.20) g(a;) = for X > 1, q^q, qeL°°[0, 1]. 
Let e (0, 1) and define the space of functions: 

(10.21) L(a;o) :=L°°(a;o,l) X L°°(A,J, 
where 

(10.22) A,„ ■.^{x,y:x(i<x<l, < \y\ < x}. 
Denote 

(10.23) ||u|| := esssup |(7(a;)| + sup \K{x,y)\. 

xq<x<1 x,y^Axo 

Let 

(10.24) \\h\\<R. 
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Theorem 10.1 (||2|). Let (10.17), (10.18), (10.20), and (10.24) hold, and choose 
any R > R. Then process (10.19) converges in L{xq) at the rate of a geometrical 
progression for any G (1 — /i, 1), where 



H := mm 



8{R-R) 2 
5R 



One has 
(10.25) 

// one starts with the data 
(10.26) 



lim Un 

n — ^CJO 



q{x) 
K{x,y) 



,x,y e A^„. 



{K{xo,y), K^{xo,y)}o<iyi<^g, 

replaces in (10.19) h by Hq :— {g°), where fo and go are calculated by formulas 
(1.17) and (1.18) in which the first argument in K{\,y),x = \, is replaced by 
X — xq, then the iterative process (10.19) with the new h ~ Hq, in the new space 



with 



L{xi) ■.^L°°{xi,xo) X L°°(A,J, 



{x,y : xi < X < xo, 0<\y\<x}, 



converges to (^K{x^y)) -^(^i)- 

In finite number of steps one can uniquely reconstruct q{x) on [0,1] from the data 
(10.12) using (10.19). 

Proof. First, we prove convergence of the process (10.19) in L{xo). 

The proof makes it clear that this process will converge in L{xi) and that in final 
number of steps one recovers q{x) uniquely on [0,1]. Let B{R) :— {U : \\U\\ < R, 
U e Lixo)}, R> R. 

Let us start with 



Lemma 10.1. The map U G W{U) + h maps B{R) into itself and is a contraction 
on B{R) if Xq G (1 - 1), min (^^^r^, 

Proof of the lemma. Let U = ( ) ; V = {k^) One has: 

2 Si {\qi " 921 \Ki\ + \q2\ \Ki ~ K2\) ds 
Ud. (l'Zi-'72||i^i| + |g2||ifi-i^2|) dsdt 



\\WiU)-WiV)\\ < 
(10.27) <\\U-V\\ 



2(l-a;o)i?+|(l-a;o)2 



< \\U^V\\il-xo)-R. 



Here we have used the estimate (1 — a;o)^ < 1 — xo and the assumption \\U\\ < R, 

< R. 
If 

2 



(10.28) 



l-Xo < 



5R 



then is a contraction on B(R). 

Let us check that the map T{U) = W{U) + h maps B{R) into itself if 1 — < pi. 
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o2 — 

Using the inequality ab < ^ if a + 6 = i?, a, 6 > 0, one gets: 

\\W{U) + h\\ < + 



(10.29) 



<T;i^-xo)^+R<R if l-a;o< ^ ' 



Thus if 

(10.30) fi = min 



2' ' 4 5i?2 
8{R-R) 2 \ 



5i?2 5R 



then the map U W{U) + /i is a contraction on B{R) in the space L{xq). Lemma 
10.1 is proved. □ 

From Lemma 10.1 it follows that process (10.19) converges at the rate of geomet- 
rical progression with common ratio (10.30). The solution to (10.14) is therefore 
unique in L{xo). 

Since for the data h which comes from a potential q S L°°{0,1) the vector 
(yK{x\)^ solves (10.14) in L{xo), it follows that this vector satisfies (10.25). Thus, 
process (10.19) allows one to reconstruct q{x) on the interval from data (10.12), 
xq = 1 — fi, where is defined in (10.30). 

If q{x) and K{x, y) are found on the interval (xq, 1), then K(xa^ y) and Kx{xo, y) 
can be calculated for < |y| < xq. Now one can repeat the argument for the interval 
(xi, xq), xq — xi < fi, and in finite number of the steps recover q{x) on the whole 
interval [0,1]. 

Note that one can use a fixed fi if one chooses R so that (10.24) holds for h defined 
by (10.16) and (10.17) with any x G L°^ [0, 1]. Such R does exist if g e i°°[0, 1]. 
Theorem 10.1 is proved. □ 

Remark 10.1. Other inverse problems have been reduced to the overdetermined 
Cauchy problem studied in this section (see p2[ , [ p2| , j. The idea of this reduc- 
tion was used in p2[ for a numerical solution of some inverse problems. 

11. Representation of I-function 

The I{k) function (2.1) equals to the Weyl function (2.3). Our aim in this section 
is to derive the following formula ([|9|): 

J • poo 

(11.1) J(fc) --ik + Y, , '""^ +a, a := a{t)e'''* dt, 

j^Q k — ikj Jq 

where ko := 0, rj = const > 0, 1 < j < J, rp > iff /(O) ~ 0, a{t) — a{t) is a 
real-valued function, 

(11.2) a{t) e L\R+) if /(O) 7^ and q{x) G ii.i(M+), 

(11.3) a{t) e L\R+) if /(O) = and g e ii,3(R+)- 

We will discuss the inverse problem of finding q{x) given I(k) Vfc > 0. Unique- 
ness of the solution to this problem is proved in Theorem 2.1. Here we discuss a 
reconstruction algorithm and give examples. Formula (11.1) appeared in |1S||. 
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Using (2.17) and (2.1) one gets 
(11.4) 

ik - A(0) + A^{y)e^^y dy 

Hk) = , /r^\,] ,' A{y):=A{0,y), A,{y) := A40,y). 

1 + Jo e^^A{y} dy 



The function 



k 1 r k '^^j 



f A- -A- 

(11.5) m = 1^1 Aiy)e^^y dy = /o(^)^ R , + ^ 

where fo{k) is analytic in C+, /o(oo) = 1 in C+, that is, 

(11.6) /o(fc) ^ 1 as ^ oo, and /o(fc) ^0, Vfc e C+ := {k : Imk > 0}. 
Let us prove 

Lemma 11.1. ///(O) 7^ and q e i/«en 

(11.7) /o(/c) = 1 + / 6o(t)e*'=* := 1 + feo, 60 e W^^'i(M+). 
Here W^^'^(IR+) is the Sobolev space of functions with the finite norm 



Polki.i := / (|6o(t)| + |6'o(i)|)dt<oo. 

JO 

Proof. It is sufficient to prove that, for any 1 < J < J, the function 
(11-8) ^^^J{k) = l + g,{t)e^'''dt, g,&W''\^+). 

Since ^4^^ = 1 + and since A{y) e W^''^{R+) provided that q e ii,i(R+) 

(see (2.18), (2.19)), it is sufficient to check that 

(11.9) = g{t)e"^' dt, gGW''\m+). 

One has f{ikj) ~ 0, thus 

Z /v J K 1 Kj J Q Ki % Kj 

poc f-y /.oc 

(11.10) =/ A{y)e-''^H e^^''-'^'^' ds = e''''hj{s)ds 

Jo Jo Jo 

where 

POO poo 

(11.11) hj{s) := i / A{y)e-''^^y-'^ dy = i A{t + s)e-^i* dt 

J 8 Jo 

Prom (ll.ll)one obtains (11.9) since A{y) € W'^'^{R+). 

Lemma 11.1 is proved. □ 

Lemma 11.2. ///(O) = Q andq& Li,2(]R+), then (11.7) holds. 
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Proof. The proof goes as above with one difference : if /(O) = then fep = is 
present in formula (11.1) and in formulas (11.10) and (11.11) with fco = one has 



(11.12) ho{s) = i I A{t + s)dt. 

Thus, using (2.18), one gets 







"'0 "'0 



OO />00 /'OO 



\ho{s)\ds<c / ds dt \q{u)\du 



t + s 



OO />00 /"OO /"OO /"OO 



(11.13) =2c/ ds dv \q{u)\du<2c ds \q{u)\udu 

Jo Jv Jo J 

poo 

= 4c / u'^ \q{u)\ du < OO if gr e Li,2(M+), 



Jo 

where c > is a constant. Similarly one checks that /iq(s) G _L^(M-|_) ii q E 
Li,2(K+). 

Lemma 11.2 is proved. □ 

Lemma 11.3. Formula (11.1) holds. 
Proof. Write 

k+i yfJ k+ikj 

, ^ 1 ^ k lij^l k-ikj 

^ ■ ^ f{k) folk) 

Clearly 

^''■''^ ^nl^-l + Er^' ^o:=0, k,>0. 

j=l J j=0 ■> 

By the Wiener-Levy theorem §17], one has 
1 1"°° 

(11.16) — - = 1 + / b{t)e'''' dt, h{t) € W^-\^+). 
fo[k) Jo 

Actually, the Wiener-Levy theorem yields b{t) G L-^(M+). 

However, since bo & W^'^{'M.+ ), one can prove that b{t) G W^'^(R+). 

Indeed, b and bo are related by the equation: 

(11.17) (l + foo)(l+^) = 1, VfcGM, 
which implies 

(11.18) b^ -bo -bob, 
or 

(11.19) b{t) = -boit) - [ bo{t- s)b{s)ds :^ -bo-bo*b, 

Jo 

where * is the convolution operation. 

Since G ii(R+) and b G iHl*+) the convolution 6o * & G LHM+). So, 
differentiating (11.19) one sees that b' G L^(]R+), as claimed. 
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Prom (11.16), (11.15) and (11.4) one gets: 
(11.20) 

m = {ik - A{0) + A,){1 + 6) I 1 + TTT 







j=0 ^ 



where c is a constant defined in (11.24) below, the constants aj are defined in (11.25) 
and the function a is defined in (11.26). We will prove that c = (see (11.28)). 
To derive (11.20), wc have used the formula: 



(11.21) ikb = ik 

and made the following transformations: 

(11.22) 



ik Jo 



h'{t)dt 



-6(0) - 6', 



I{k) = ik- A(0) - 6(0) - 6' + - A(0)6 + AibY, 



^ c,[^(o)+6(o)] ' m-mi) 
^ k-ikj • ■• 



j=0 



k ikq 

j=0 ^ 



^ k- iki '' 



where 
(11.23) 



g{k) := -b' + - A{0)b + A^b. 



Comparing (11.22) and (11.20) one concludes that 



(11.24) 
(11.25) 
(11.26) 



c := 



-A(0)-6(0)+z^c„ 
-Cj[kj+A{Q) + b{Q)-g{ikj)], 



k - itj 
j=o J 



a. 



To complete the proof of Lemma 11.3 one has to prove that c = 0, where c is 
defined in (11.24). This is easily seen from the asymptotics of I{k) as k oo. 
Namely, one has, as in (11.21): 



(11.27) 



ik 



ik 



Prom (11.27) and (11.4) it follows that 

^(0) 



I{k) = {ik-A{0) + Ai) 



ik 



(11.28) ={ik-A{0)+Ai) (^1 + ^+0^^ 

Prom (11.28) and (11.20) it follows that c = 0. 
Lemma 11.3 is proved. 



= ik + o{l), fc — > +00. 



□ 



Lemma 11.4. One has aj = irj, rj > 0, 1 < j < J , 
ro>Oi//(0) = 0. 



ro = z//(0) 0, and 
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Proof. One has 



(11.29) flj ^ Res I{k) 

From (2.7) and (11.29) one gets: 



f{ikj 



c 



ZjOiXij J 3 



If j = 0, then 



(11.31) ao = Res/(fc) 

k—0 



/'(o,o) 



/(O) 

Here by Resfe^o we mean the right-hand side of (11.31) since I{k) is, in general, 
not analytic in a disc centered at fc = 0, it is analytic in C4. and, in general, cannot 
be continued analytically into C- . 

Let us assume q{x) £ Li,2(IR-|_). In this case f{k) is continuously differentiable 
in C+. 

From the Wronskian formula 

(11.32) mm-k)-f'io-k)m^^^ 

k 

taking fc — > 0, one gets 

(11.33) /'(0,0)/(0) = -z. 

Therefore if g G ii,2(IR+) and /(O) = 0, then /(O) ^ and /'(0,0) 7^ 0. One 
can prove ^ pp.188-190], that if g G Li,i(K+), then is bounded as fc 0, 
fc G C+. 

From (11.31) and (11.33) it follows that 



(11.34) aj = 7T = zro, ro := 

From (2.17) one gets: 

(11.35) /(0)=«/ A{y)ydy. 



Since A{y) is a real-valued function if q{x) is real-valued (this follows from the 
integral equation (5.62), formula (11.35) shows that 

(11.36) [/(0)]'<0, 
and (11.34) implies 

(11.37) ro > 0. 

Lemma 11.4 is proved. □ 

One may be interested in the properties of function a{t) in (11.1). These can be 
obtained from (11.26), (11.16) and (11.7) as in the proof of Lemmas 11.1 and 11.2. 
In particular (11.2) and (11.3) can be obtained. 

Note that even if q(x) ^ is compactly supported, one cannot claim that a(t) is 
compactly supported. 

This can be proved as follows. 
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Assume for simplicity that J — and /(O) 7^ 0. Then if a{t) is compactly 
supported then I{k) is an entire function of exponential type. It is proved in 
p. 278] that if q{x) ^ is compactly supported, q G L^(]R+), then f{k) has infinitely 
many zeros in C. The function /'(O, z) ii f{z) — 0. Indeed, if f{z) = and 
f'{0,z) = then f{x,z) = by the uniqueness of the solution of the Cauchy 
problem for equation (1.1) with k ^ z. Since f{x,z) ^ (see (1.3)), one has a 
contradiction, which proves that /'(O, z) 7^ if f{z) — 0. Thus I{k) cannot be an 
entire function if q{x) ^ 0, q{x) G L^(R+) and q{x) is compactly supported. 

Let us consider the following question: 

What are the potentials for which a{t) = in (11.1)? 

In other words, suppose 

,7 

(11.38) I{k)^ik + Y^ 



7=0 ^ 

find q{x) corresponding to /-function (11.38), and describe the decay properties of 
q{x) as a; — > +00. 

We now show two ways of doing this. 

By definition 

(11.39) f{0,k)=I{k)f{k), f{0,-k) = I{-k)f{-k), keR. 

Using (11.39) and (2.23) one gets 

[I{k)-Ii-k)]fik)fi-k) = 2ik, 



(11.40) /(A:)/(-fc) = VfceM. 

By (2.5), (2.6) and (11.30) one can write (see (l9|) the spectral function corre- 
sponding the /-function (11.38) (VA = k): 

(11.41) *(A) JT""'' , . 

where 5{X) is the delta-function. 

Knowing dp{X) one can recover q{x) algorithmically by the scheme (5.26). 
Consider an example. Suppose /(O) ^0, J = 1, 

/-,, .r.s ^"f"! -1 iri{k + iki) . /, rik \ rik 

(11.42) I{k)^ik+- —^ik-\ -p, T2^^^\k+ ' 



k - iki fc2 ^ y fc2 + fcf / fc2 + fcf 

Then (11.41) yields: 



(11.43) dp{X) 
Thus (5.27) yields: 



2kiriS{X + k()dX, A<0 



1 /"^ nVA sin\/Axsin%/Aj/ sh{kix) sh{kiy) 

(11.44) L(x,y) ^ - I dX^—i^ — — h2fciri- 



TT 







A + /c? VA \/A ki ki 
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and, setting X = and taking for simplicity 2fciri — 1, one finds: 

dkk^ sin(fca:) sm.{ky) 

° d/c sin(fca:) sin(A;?/) 

dk[cos k{x — y) — cos fc(a; + y)] 



(11.45) 



Lo[x,y) := 

TT 

= ^ 
TT 

~ 2fci 



fc2 + fc2 



where the known formula was used: 



(11.46) 

Thus 

(11.47) 



cos kx 

k"^ + a? 



dk 



-a\x\ 



2a 



a > 0, 



fci > 0, 



L{x,y) 



2ki 



-ki\x-y\ 



-ki(x+y) 



sh{kix) sh{kiy) 



ki 



Equation (5.30) with kernel (11.47) is not an integral equation with degenerate 
kernel: 

"g-fci|t-j/| _g-fci(t+t/) 



(11.48) K{x,y)+ / K{x,t) 



sh{kit) shikiy) 



2fci/ri fci 

^-k^\x-y\ _ g-fci(x+s,) sh{kix) Sh{kiy) 



dt 



fci 



fci 



2fci/ri 

This equation can be solved analytically jS^ , but the solution requires space to 
present. Therefore we do not give the theory developed in |Q but give another 
approach to a study of the properties of q{x) given I{k) of the form (11.42). This 
approach is based on the theory of the Riemann problem 

Equations (11.40) and (11.42) imply 

e + ki 



(11.49) 
The function 
(11.50) 

Write (11.49) as 

Thus 
(11.51) 



f{k)f{^k) 



fc2 



fcj + ri. 



/o(fc)-/(fc)^^0 in 



/o(fc)^^/o(-fc) 



fc — iki 

k + ifci 



e + ki 



/o(fc) 



k + iki 

e + ki 



k — iki fc^ + 



h 



h{k) 



1 



fc2 + ^2 

The function fo{—k) 7^ in C_, /o(oo) = 1 in C_, so h 
Consider (11.51) as a Riemann problem. One has 



M-k)- 

f^(^_i;^ is analytic in C- 



(11.52) 



inda 



e + ki 

fc2 + vl 



1 

27ri 



din 



kl 



fc2 



0. 
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Therefore (see problem (11.51) is uniquely solvable. Its solution is: 

as one can check. 
Thus, by (11.50), 



(11.54) /(fc) 



k — iki 
k + ivi 



The corresponding S'-matrix is: 
, . f{-k) _ {J^ + zh)ik + u^ 

^ ^^^^ ^ ^" W(fc-^fci)(fc-*-i) 

Thus 

1 f°° 

(11.56) Fs{x):^— [1- S{k)]e'''''dk^O(e-^'=') for a; > 0, 

Frf(a;) = si e-'=^^ 

and 

(11.57) F{x) = Fs{x) + Fd{x) = O (e^'^^-) . 
Equation (5.50) implies A{x,x) = O (e'^''^^), so 

(11.58) g(x) = O (e-2fei^) , x ^ +oo. 



Thus, if /(O) 7^ and a{t) = then q{x) decays exponentially at the rate 

termined by the number k^, k^ = min kj. 

i<j<J ' 

If /(O) = 0, J = 0, and a{t) = 0, then 



(11.59) /(fc) = ifc + ^, 

A; 

(11.60) f{k)fi~k) - -2^, ro > 0. 
Let /o(fc) = it+lKLtl^ Then equation (11.60) implies: 

(11.61) Mk)M-k) - := ^0, 

and fo{k) ^ in C+. 

Thus, since inds. ^2 — Oi fo{k) is uniquely determined by the Riemann prob- 
lem (11.61). 



One has: 



and 

(11.62) m-T^^ ^(^) = 47#-|^' 



27rJ_ooV k-ivoj 27r J _^ k - ivq 
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and Fd{x) = 0. 
So one gets: 

(11.63) F{x) = Fsix) = 2uoe-'""', x > 0. 
Equation (5.50) yields: 

POO 

(11.64) A{x,y) + 2vQ I A{x,t)e-''°^^+y">dt = -2voe-''°^''+y\ y>x>0. 

J X 

Solving (11.64) yields: 

(11.65) A{x^y) = -2z.oe--(^+^)^^^^. 

The corresponding potential (5.51) is 

(11.66) q{x) = O , X ^ 00. 

If q{x) — 0{e^'^^), k > 0, then a(t) in (11.1) decays exponentially. Indeed, in 
this case b'(t), Ai{y), h{t), Ai * b decay cxpcncntially, so, by (11.23), g{t) decays 
exponentially, and, by (11.26), the function ^'"'^IZ^^^'^'' ^ '■= h with h{t) decaying 
exponentially. We leave the details to the reader. 

12. Algorithms for finding q{x) from I{k) 

One algorithm, discussed in section 11, is based on finding the spectral function 
p{\) from I{k) by formula (11.41) and then finding q{x) by the method (5.26). 

The second algorithm is based on finding the scattering data (2.10) and then 
finding q{x) by the method (5.49). 

In both cases one has to find kj, 1 < j < J, and the number J. In the second 
method one has to find f{k) and Sj also, and S{k) = ^j^- 

If kj and f{k) are found then sj can be found from (2.12). Indeed, by (11.1) 

(12.1) irj := Res I{k) = Z^i^llM. 
Prom (12.1) and (2.12) one finds 

(12.2) sj= ' ~ 



If kj arc found, then one can find f{k) from I{k) as follows. Since f'{0,k) = 
f{k)I{k), equation (2.23) implies equation (11.40): 

(12.3) = _L_. 

Define 

(12.4) -W:=n|^ if /(0)<oo, 
and 

(12.5) w{k):=-^r\'^^^ if 7(0) =oo. 

^ ^ ^ ^ k + i k + iki ^ ^ 

One has 1(0) < oo if /(O) and /(O) = oo if /(O) = 0. Note that if g G ii,2(R+) 
and /(O) = then /'(0, 0) 7^ and /(O) ^ 0. 
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Define 

(12.6) h{k) := ffl. 

w(k) 

Then h{k) is analytic in C_|_, h{k) ^ in C+, and h{oo) = 1 in C+, while h{—k) 
has similar properties in C_. Denote ^^j'^^^ := h^{k). This function is analytic in 

C_, h^{k) 7^ in C_ and h^{oo) = 1 in C_. Denote h{k) h+{k). 
Write (12.3) as the Riemann problem: 

(12.7) h+{k) ^ g{k)h^{k), 
where 

(12.8) SC^^-J^) if ^W<-' 
and 

(12.9) ^^(^)-7;^^ if m-^- 

We claim that the function g{k) is positive for all A: > 0, bounded in a neighborhood 
of fc = and has a finite limit at fc = even if /(O) = 0. Only the case /(O) = 
requires a comment. If /(O) = 0, then /'(0,0) = 0, /'(0,0) ^ 0, /(O) ^ 0, and 
one can see from (12.3) that the function is bounded. Thus, the claim is 

verified. 

The Riemann problem (12.7) can be solved analytically: ln/i+(fc) ~ lnft,_(fc) = 
In g{k) and since h+{k) and h_ (fc) do not vanish in C+ and C_ respectively, In h^{k) 
and ln/i_(fc) are analytic in C+ and C_ respectively. Therefore 



(12.10) 



(12.11) h{k)^h+{k) if /mfc>0, h{k) = h^{k) if /mfc<0, 
and 

(12.12) f{k) = w{k)h{k), lmk>0. 

Finally, let us explain how to find kj and J given /(fc). 
From (11.1) it follows that 

1 r 

(12.13) — J (/(fc) - ifc)e-*'=* dk = ~Y2 ^o^^'^ - y fo^' ^ < 0- 

Taking t — s- — oo in (12.13) one can find step by step the numbers tq, fci, ri, fc2, 
r2 ■ ■ ■ , rj, kj. If /(O) < oo, then ro = 0. 

13. Remarks. 

13.1. Representation of the products of the solution to (1.1). In this sub- 
section we follow |l^. Consider equation (1.1) with q = qj, j = 1,2. The function 
u(x, y) := ipi{x, k)ip2{y, fc) where ipj, j — 1, 2, satisfy the first two conditions (1.4), 
solves the problem 

\ 1 r 92 

(13.1) a^"^^*'^'' ^""^2(2/) u{x,y), 
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(13.2) u{0, y) = 0, u^(0, y) = ^y, k), 

(13.3) u{x,0) = 0, Uy{x,0) = tpi{x,k). 
Let us write (13.1) as 

(13-4) ^ - ^ ^i^' y) = bi(a:) - q2{y)Hx, y) 



\ dx^ dy'^ ^ 

and use the known D'Alembert's formula to solve (13.3)-(13.4): 

If 1 r+y 

(13.5) u{x,y) = - [qi{s) - q2{t)]u{s,t)dsdt+ - ipi{s)ds, 

J D^y ^ J x—y 

where D^y is the triangle Q<t<y,x — y + t<s<x + y — t. 

Function (13.5) satifies (13.3) and (13.1). Equation (13.5) is uniquely solvable 
by iterations: 

°° I rx+y 

(13.6) u{x,y) = ^Um{x,y), uo{x,y) := - ipi{s,x)ds, 

(13.7) Um+i{x,y) = ^ f [qi{s) - q2{t)]um{s,t)dsdt. 



Note that 



I i-x+y 

(13.8) Um{x,y) = - Wm{x,y,s)tpi{s)ds. 

^ Jx-y 

If TO = this is clear from (13.6). If it is true for some to > 0, then it is true for 
TO+ 1: 



1 rV rx-ty-z ^ ns+t 

Um+i{x,y) = - dt I ds[qi{s) - q2{t)]- I Wm{s,t,a)(pi{a) da 

^ Jo Jx-y+t ^ Js-t 

I rV rx+y ^ 

(13.9) =0/ / da<fii{a)wmix,y,t,a) 

^ Jo Jx-y 
I rx+y 

= -/ da(fii{a)wm+i{x,y,a), 

^ Jx-y 

where Wm and Wm+i are some functions. 

Thus, by induction, one gets (13.8) for all to, and (13.6) implies 

I rx+y 

(13.10) u{x,y)=^- w{x,y,s)(pi{s)ds, 

Jx-y 

where 

00 

(13.11) w{x,y,s) := ^Wm{x,y,s). 

m=0 
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To satisfy (13.2) one has to satisfy the equations: 

0= f w{0,y,s)(pi{s,k)ds, 

J-y 

(13.12) V2iy,k) = ^[w{0,y,y)ipi{y) - w{0,y,-y)ipi{-y)] 

1 f" 

+ 2 J '^x{0,y,s)(pi{s)ds. 

Formula (13.10) yields 

(13.13) (fii{x,k)(fi2{y,k) = - w{x,y,s)(fii{s,k)ds. 

^ Jx-y 

If X = y, then 

1 /■^'^ 

(13.14) (fii{x,k)ip2{x,k) = - / 'w{x,x,s)(pi{s,k)ds. 

2 Jo 

Therefore, if 

r 

/ h{x)ipi{x,k)(p2{x,k) dx = Vfc > 0, 
Jo 

then 

pa p2x 

0= / h{x) I w{x,x, s)ipi{s,k) ds dx 
Jo Jo 

p2a pa 

= I dsipi{s,k) / dxh{x)w{x,x,s) VA; > 0. 
Jo Ji 

Since the set {(pi{s, fc)}vfe>o is complete in 1/^(0, 2a), it follows that = J" dxh{x)w{x, x, s) 
for all s e [0, 2a] . Differentiate with respect to s and get 

(13.15) ^ (I' f (i) ~ dsh{x)ws{x,x,s) =Q. 

From Voltcrra equation (13.15) it follows /i(.x) = if the kernel u>s (a;, a;, s)u>^^ (f'f'^) 
t{x, s) is summable. From the definition (13.11) of w it follows that if \<l{x)\ dx < 
oo V6 > 0, then Ws{x,y,s) is summable. The function w{x,y,s) has m summable 
derivatives with respect to x, y and s if q{x) has m—1 summable derivatives. Thus 
one can derive from (13.15) that h{x) = if w (|, |, s) > for all s S [0, 2a]. 

If the boundary conditions at a; = are different, for example, ip'-{Q,k) — 
hoiPj{0,k) = 0, j = 1,2, then conditions 

(13.16) Ux — hou\_^_^ = 0, Uy — hou\^_^ = 0, ho = const > 
replace the first conditions (13.2) and (13.3). One can normalize <^j(a;, k) by setting 

(13.17) (pj{0,k) = l. 
Then 

(13.18) 'f'j{0,k) = ho, 

(13.19) u{0,y) = ip2{y,k), u{x,0) = ifii{x,k), 

(13.20) Ux{0,k) = hoip2{y,k), Uy{x,0) = ho(pi{x,k), 
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and (13.5) is replaced by 



'^{^^y) = 7: / [qi{s) ~ q2{t)]u{s,t)dsdt 
^ JD^y 

1 r+y 1 

(13.21) hoLpi{s,k)ds + -[ipi{x + y,k) + (fi{x - y,k)]. 
Note that 

1 Id /■^+^ 1 f^^y 

(13.22) -^[Viix + y) + ^lix - y)] ^ -— J ipi{s) ds ^ - J (psis)ds. 

Equation (13.21) is uniquely solvable by iterations, as above, and its solution is 
given by the first formula (13.6) with 

ho /"^^^ 1 

(13.23) uo{x,y)^— ipi{s,k) ds + - i ipis{s,k)ds. 

J x — y J x — y 

The rest of the argument is as above: one proves existence and uniqueness of 
the solution to equation (13.21) and the analog of formula (13.10): 



1 



x+y 



(13.24) u{x,y)^-l w{x,y,s)^{s)ds, ^ hofi{s,k) + ipisis,k). 

^ J x—y 

Thus 

(13.25) u{x,x) = (pi{x,k)ip2{x,k) ^ I t[x, s)[hoLpi{s,k) + ipis{s,k)\ds, 

Jo 

where 

(13.26) t{x,s) :— -w{x,x,s), 

and t(x, s) is summable. 

Thus, as before, completeness of the set of products {ipi(x, k)ip2{x, k)}yk>o can 
be studied. 

13.2. Characterization of Weyl's solutions. The standard definition of Weyl's 
solution to (1.1) is given by (2.2). 
In [|l^ it is proved that 

(13.27) T4^(2:,A:) = e''=^(l + o(l)) as \k\ ^ oo, \x\ < b, fc^ e A, 
where A := {A : \Im\\ > e, dist{\,S) > e}, 

(13.28) 5 := MU [i7_,i7+], 7± := inf ± [u'^ + q\u\^] dx . 

u(0)=0 

The relation (13.27) gives a definition of the Weyl solution by its behavior on 
compact sets in the i-space as oo, as opposite to (2.2), where k is fixed, and 

X — > cxD. For multidimensional Schrodinger equation similar definition was proposed 
in 0, p.356, problem 8]. 

We want to derive (13.27) for potentials in Li_i(M_(.) and for k > 0, k ^ +oo. 
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The idea is simple. For any q = q G Ll^^{M.-^-), one can construct ip{x,k) and 
^{x, k), the solutions to (1.1) and (1.4), for any |a;| < b, where 6 > is an arbitrary 
large fixed number, by solving the Volterra equations 

/ ,N sin(fca:) sinlkix — y)] , , , , 

(13.29) ifix, k) = — ^ + ^ ^'^ q{yMy, k) dy 

(13.30) ^(X, k) = COs(fcx) + r 51!lM^Jg(y)^(y, fc) dy. 



k 

One can also write an equation for the Weyl solution W: 
(13.31) W{x, k) = cos(fca;) + m(fc)^^4^ + f ''""^^^f ~ ^\ (y)W^(2/, k) dy. 



This equation is uniquely solvable by iterations for < h. 
It is known that 

(13.32) m(/c) = iA: + o(l), ^ oo, Imk > e\Rek\, e > 0. 

For q E Li,i(M+) the above formula holds when fc > 0, fc ^ +oo. From (13.31) and 

(13.32) one gets, assuming A; > 0, 

(13.33) W{x, k) = (^1 + O ) + fc) dy. 

Solving (13.33) by iterations yields (13.27) for fc > 0, fc ^ +oo. For q e Li4(M+) 
the Weyl solution is the Jost solution. Therefore the above result for k > 0,k ^ +oo 
is just the standard asymptotics for the Jost solution. It would be of interest to 
generalize the above approach to the case of complex k in the region (13.27). □ 
One can look for an asymptotic representation of the solution to (1.1) for large 
|fc|, Imk > e\Rek\, e > 0, of the following form: 

(13.34) u{x, k) = e*'=^^+i'o" 



where 

(13.35) a' + 2ik(7 + (t'^ - qix) =0, cr = ^ + o | ^ ) , |/c| -> oo. 

2ik \k J 

From (13.34) one finds, assuming q{x) continuous at x = 0, 

If q{x) has n derivatives, more terms of the asymptotics can be written (see [ p^ 
p.55]). 

13.3. Representation of the Weyl function via the Green function. The 

Green fvmction of the Dirichlet operator Lq = —j^+q{x) in i^(]R+) can be written 
as: 

(13.37) G{x,y,z)^ip{y,^z)W{x,^z), x>y 

where ip{x,k), k := y/z, solves (1.1) and satisfies the first two conditions (1.4), and 
W{x, y/z) is the Weyl solution (2.2), which satisfies the conditions: 

(13.38) W{0, yi) = 1, W'{0, VI) = m{yG). 
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From (1.4), (13.37) and (13.38) it follows that: 

If g e ii,i(K+) then W{x,k) = ^j^, where f{x,k) is the Jost solution (1.3), 
k e C+. Note that (13.17) and (13.38) imply: 

(13.40) hmG{x,y,X + iO) = Jmm(7AT^) 

TT TT 

and 

(13.41) G{x, y, z) = r ^^iM dp{t), e{x, y, t) = cp{x, Vi^y, Vt). 
Thus 

(13.42) -Im G{x, y, t + iQ)dt = e{x, y, t)dp{t). 

TT 

From (13.42) and (13.40) one gets, assuming p(— oo) = 0, 

1 ft ^ 

(13.43) p{t) = - Imm{VX + iO) d\. 

7^ J-oo 

If A < then Imm (\/A + iO) = except at the points A = — at which 



f{ikj) = 0, so that m{y/ —k"^ + iO) = oo. Thus, if t and a are continuity points of 
p{t), then 



1 



(13.44) p{t)-p{a) = - Imm{VxTiO)dX, a > 0. 

J a 

Let us recall the Stieltjes inversion formula: 

li z = a + IT, T > 0, pit) is a function of bounded variation on M, 

(13.45, -w-r?^. 

y_<x> t-Z 

and if a and b are continuity points of p{t), then 
1 

(13.46) - / Innp{X + iQ)dX = p{b) - p{a). 

J a 

Therefore (13.44) implies 

(13.47) miV-z) - r 

J — CO Z 

The spectral function dp{t) does not have a bounded variation globally, on the 
whole real axis, and integral (13.47) diverges in the classical sense. We want to 
reduce it to a convergent integral by subtracting the classically divergent part of it. 

If q{x) = 0, then p := po{t) for t < 0, m(-\/A) = iVX, and formula (13.44) with 
a = yields 

2Ai 

(13.48) po{X) = —. 
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If q{x) = then G{x, y, A) = Sill^Mi e^^^, y < x, so (13.39) yields m{VX) = 
iy/X. Formula (13.47) yields formally 

(13.49) zVA^if^. 

TT Jo t- X 

This integral diverges from the classical point of view. Let us interpret (13.49) as 
follows. Let ImX > 0. Differentiate (13.49) formally and get 

(13.50) J^ = ir^^, 7mA >0. 

This is an identity, so (13.49) can be interpreted as an integral from to A of 
(13.50). The integral t' 2 dt which one obtains in the process of integration, 
is interpreted as zero, as an integral of a hyperfunction or Hadamard finite part 
integral. 

Subtract from (13.47) the divergent part (13.49) and get: 
(13.47') m{V-z)-iV~z= r 

where 

{ Vxdx \ > n 
U, A <C U. 

Integral (13.47') converges in the classical sense ii q E Li.i(K+). Indeed, by (2.5) 
and (13.47") one has da{t) = ^ijJl^^ - ^)dt. By (5.65) one has /(V<) = 1 + 
O(^) as t +00. Thus da{t) = 0{^)dt as t +00. Therefore integral (13.47') 
converges in the classical sense, absolutely, if 7m^; ^ 0, otherwise it converges in 
the sense of the Cauchy principal value. 
Let us write (11.1) as 

(13.51) 

J 



m{k) — ik = J 



00 



dt, H{t) 



1, t>0, 
0, t<0. 



From (13.51) and (13.47') one gets 

(13.52) r = r e"=*a(t) dt, X = k^ + iO, 
J —00 ^ ^ J —00 

where 

.7 

(13.53) a{t) := - XI ^ie'=-'*77(-t) + a{t)H{t). 

j=o 

Taking the inverse Fourier transform of (13.52) one can find a(t) in terms of 
a{s). If A: > then k = Vk^ + iO and if fc < then k = Vfc^ - iO- Thus: 

1 r ...-^kt r dais) 



ait) 



/CO pC 
dke-'''' / 
-CSO J — ( 



1 fX r-oc -ikt f 



kt pO ^—ikt 

dk 
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Let us calculate the interior integral in the right-hand side of the above formula. 
One has to consider two cases: s > and s < 0. Assume first that s > 0. Then 



JO 



-ikt 



dk 



k'^ + iO-s 



-ikt 



dk 



k"^ — iO — s 



ikt 



dk 



f 

Jo 



e-'''*S{k'^-s)dk- 



-ikt 



(13.55) 
where 

If s < 0, then 
(13.56) 

If s > 0, then 
(13.57) 

From (13.54)-(13.57) one gets 
(13.58) 



5(A;2 - s)dk 



■-r 

J — C 



— -=^\VLU~st)H{s) + J, 
^/s 

-'''*dk 



J = 



k'^ — s ' 



-\/HI*I 



J 



sin(|f|^/s). 



f°° , , , sm(tJs) , If 
a{t) = da{s) ^ ' H{t) - - y 



Formula (13.58) agrees with (13.53): the second integral in (13.58) for t > is 
an L'^(R+) function, while for t < it reduces to the sum in (13.53) because 
da{s) = dp{s) for .s < 0, dp{s) for s < is given by formula (2.5) and the relation 
between cj and rj is given by formula (2.7). 
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